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Nonsteady Supersonic Flow About 
of Revolution 


Pointed Bodies 


WILLIAM H. DORRANCE* 
Unwersity of Michigan, Willow Run Research Center 


(I) SUMMARY 


The general velocity potential for a slender pointed body of 
revolution oscillating periodically in a supersonic uniform stream 
is presented. The solution is shown to reduce in the limit to the 
sum of the well-known stationary source-sink and doublet solu- 
tions when the frequency of oscillation approaches zero. Sta- 
bility derivatives for low-frequency oscillations are determined by 
making use of the potential expression obtained by expanding the 
oscillating velocity potential in powers of reduced frequency, 
retaining only the first-order terms in reduced frequency. In 
particular, the following four cases of steady and nonsteady 
(1) harmonic pitching about a 
(3) steady 


supersonic motion are treated: 
point x; (2) harmonic normal oscillations in pitch; 
pitching about a point x»; and (4) steady angle of attack 


tions for the slender body theory normal force and pitching mo 


Equa 


ment stability derivatives are also determined. 

Two types of bodies of revolution are examined to demon 
strate how the aerodynamic stability derivates may be varied 
by changing the contour of the body. ‘“‘Flaring’’ the aft end of 
the body is shown to increase the lift, as well to increase the 
aerodynamic damping moment of a finless body. 


(II) SYMBOLS AND NOMENCLATURE 


local velocity of sound 


normal force coefficient 


pitching moment about x» (positive in stall) 

OC),/O( el / Us), OCm/O( gl. / Uo) | dimensionless non- 
steady pitching stability derivatives | 

OC, /AwL/Us?), OCm/A(wL/ Uo?) 
normal acceleration stability derivatives] 

OC,/A(qL/Uo), OCn/%XgL/ Uo) 
steady pitching stability derivatives] 

OC, /Oa, OCm/Oa 
cients] 

—2[(¢2,/U0) + (¢2,/Uc?)] [incremental 


{dimensionless 
[dimensionless 


[steady aerodynamic coeffi- 


lifting 
pressure coefficient | 

exp | —sx} f(x) dx [Laplace transtorm of f(x)] 
0 


modified with 


K((z = Bessel function of second kind 
modulus z 


Received November 20, 1950. Revised and received March 
2, 1951. 


P —_— 
Aerodynamics Group Leader. 


(M/B) (wr/ Ue) |reduced frequency] 

length of body of revolution 

free-stream Mach Number 

local static pressure 

free-stream dynamic pressure 

steady pitching velocity (Fig. 2c) 

locus of surface of body of revolution 

reference area 

cross-sectional area of body of revolution at sta- 
tion x 

time, sec. 

free-stream speed 

local velocity 

upwash velocity (Fig. 1 

center of rotation of pitching motion 

cylindrical coordinates 

body stability axes through x 

angle of attack 

V/ Me — 1 compressibility parameter 

ratio of specific heats 

free stream velocity potential 

perturbation velocity potentials 

total velocity potential 

frequency of periodic motion 


(III) INTRODUCTION 


} emery HAS BEEN DIRECTED lately to the prob- 
lem of dynamic stability of low aspect ratio high- 
speed aircraft.’ * 
static and dynamic aerodynamic force and moment 
stability derivatives of low aspect ratio wings at super- 
sonic speeds has been treated by many investigators.* 

Little consideration beyond first approximations has 
been given to the corresponding problems for bodies of 
revolution.’ The dynamic stability derivatives for 
such bodies become of prime importance for finless 
missiles and missile configurations having relatively low 


The problem of determining the 


area, low aspect ratio surfaces. 

It is the purpose of this paper to present the general 
potential for the harmonic motion of a slender pointed 
body of revolution in a steady supersonic stream. Fol- 
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lowing this, the doublet potential for low-frequency har- 


monic oscillations is demonstrated. 
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tions (except as a product with the g potential velocity 
The potential g2, being odd with re- 


perturbations). 
spect to the plane 6 = 7 2, will yield the desired lifting 


The slender body 
theory doublet potentials are also determined. Ex- 
ample solutions for a body of revolution are calculated velocity perturbations. 
to illustrate the variation of the stability derivatives 
with Mach Number in a limited supersonic Mach Num- 
ber range. It is shown how the stability derivatives can 
be modified considerably by varying the contour of the 


(a) The Nonsteady Source-Sink Potential 


Let 
body of revolution. 

The solutions for steady pitching and steady angle 
of attack are also presented for comparison purposes. 


gi = —n(x, 7) exp jtw(t — (Mx/cB?)]} (5 


in Eq. (3). There results 


l w" 


FUNDAMENTAL SOLUTIONS Mr + ~~ B°nrs — 2347 
g c*; 


II 


(IV) 


The well-known linearized potential equation for ,, ; : 
. os i‘ ” 6 Further, make use of the Laplace transformation, and 
compressible nonsteady flow is, in cylindrical coordi- 


nates (Fig. 1), let 

: l l l 2M — i 

Brrr 7 aa 7 re Yoo — Yr + ott + Gri = O (1 ) nS, i} = exp = SX { 
r Cc” c J 0 


r- 


n(x, r)dx (j 


The solutions with which this paper is concerned will 


satisfy Eq. (1). _As such, these solutions are all subject n(0, r) = (On dx) (0, r) =0 

to the usual limitations of the linearized theory. The 

solution sought in this paper is appropriate to pointed bodies, operating upon Eq. (6 

with Eq. (7) yields 
yo = ¢1, COS 6 (2) 
; d’y 1d oa oh. 
where y; is the solution to the equation a ris +r je = 0 s 
ar* a7 cB 
| I 2M - : : a 
Born — Orr — ~ Ort Ont ¢r, = 0 (3) The desired solution to Eq. (8) is the modified Bessel 

, ail i function of the second kind. 


Direct substitution of Eq. (2) into Eq. (1) will verify 4 = Ko{ Br[s? + (w?/c284) |" (9 


this. Here, the perturbation potential representing 
the body is composed of Using transform pair (871.5) of reference 11 we obtain 
i aa ge =" cos [(w/cB?) (x? — B2r?)’ 
¥ gi + ¢2 (4) n(x, r) = | : ft | 
— P . ‘ ° P (x* — 6*r*)”* (10 
rhe potential ¢, is axially symmetric with respect to 
the x-axis and cannot yield “‘lifting’’ velocity perturba- x — Br>0 
The point source solution to Eq. (3) is obtained through Eqs. (5) and (10). 
exp jiw [t — (Mx cB?) |} cos [(w/cB?) (x? — B2r*) | 
ci (x? — B?r?)’ (11 
x— pF > DO 
If the sources are distributed along the x-axis at points £ according to f(&), then Eq. (11) becomes 
"%—Br tw M \ jw 
, f(é) exp 4 — a E) e cos . [iw — &)? — B'r?] 7¢ dé 
i = —exp } it | a | cp wie cB? f “| , 
J 0 (x — £)? — Br? , | = 
x — Br>é 
If w = 0, asin the steady source-sink solution, then Eq. (12) becomes the easily recognized form 
(13) 


<i ws 


aad f(&dé 
- o x 
0 I(x et é )2 ee B?r? 





ocity 


h re- 
fting 


and 


tain 
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(b) The Nonsteady ‘*‘Doublet”’ Potential 


Eqs. (2) and (12) are now employed to obtain the desired supersonic nonsteady doublet potential. 


entiation is accomplished by first substituting 


This differ- 


u = cosh! [(x — &)/6r] 


to eliminate the improper integral, which is then easily differentiated and results in Eq. (24). 


verted to the original form by 


so that the resulting doublet potential is 


» = exp}twt{ cos 6 = 
rcp a 


ra By f 
f'(&) (x — &) exp jo 


ms vy... j 
pAx = ep >= 
iwM He) al 


cB" ‘ 0 


If, in Eq. (14), » = O, the recognized steady doublet 


potential results: 


=i f° Br f'(E) (vw — Eddé | 
Sor oy (x — &)? — B'?]") (45) 


x — Br>€é | 


Ey. (14) is the desired nonsteady doublet solution to 
Eq. (1). The balance of this paper will be devoted to 
presenting particular solutions to Eq. (1) appropriate 
to: (1) harmonic pitching about a point x»; (2) har- 
monic normal oscillations in pitch; (3) steady pitching 
about a point xo; and (4) steady angle of attack. All 
of these situations can be treated by making use of 


Eq. (14) and its derivatives. 


¢: Boundary Conditions 

The boundary conditions for the steady or nonsteady 
supersonic flow about a body of revolution specify that 
the velocity component normal to the surface of the 
body at the surface must be zero. That is, if 


r = R(x) Qs x S34) (16) 


is the locus of the surface of the body in all meridian 
planes, then the boundary condition is (Fig. 1) 
(Ogve/Or), =-= Ro = —W(X, t) cos 6 (17) 

in general. Eq. (17) takes a different form for each of 
the four steady and nonsteady motions previously men- 
tioned. 

(1) Harmonic Pitching About the Point xo.—In this 
case 


a(t) = a exp | iwt { 


This may be con- 


£& = x — Br cosh u 


- ~x — BF si j eM all . fo - —— it 
/ f() exp 4 — ty £) ¢ sin ge (@ — §)? — B*r?] (dé + 


tw A mo w in 
“8? (vw — &)¢ cos a [(w — &)? — B?r?] dé 
r [(x — £)? — B?r?]’ T 
a, hei a | 
cB E) ¢ cos ca? [(x £) B?r?| dé a 
(x — £)? — Br? ] 


and from Fig. 2a and Eq. (17) 


(**) _ 
Or R—~0 


_ ; Unan + tw(x — Xo) ao4 exp | twot | cos 6 (18a) 


(2) Harmonic Normal Oscillations in Pitch.—From 


Fig. 2b and Eq. (17), 
(O¢g2/Or) rp +o = —Wo exp }iwt} cos 6 (1Nb) 


(3) Steady Pitching About the Point xo.-From Fig. 


2c and Eq. (17), 


(Oge/Or)r +9 = —aol/o cos 6 — g(x — Xo) cos@ (18e) 


The capital letters X and Z in Fig. 2c denote the body 
stability axes. The point xo is the center of gravity, 
and the steady pitching velocity is taken in the con- 
ventionally defined sense. 


(4) Steady Angle of Attack.—From Fig. 2d and Eq. 
(17), 


(Og2/Or)r +o = —Uo a cos 0 (1Nd 


This case has been treated in the literature'* '* and 
will not be dwelt upon in any great detail in this paper 
other than to compare it with the other solutions. 

Eqs. (14) and (15), together with Eq. (1S), are now 
used to specify the function f(g). The value of (O¢ge 
Or)r—o is obtained by the same limiting process as in 


references 12 and 13 so that 
y2/dr)p-o = [—fle)/R?] exp fiat} cos (19 


The function f(g) is now determined through Eqs. 
(18) and (19) for all of the cases treated in this paper. 
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(V) NorMAL ForCE AND PITCHING MOMENT 


COEFFICIENTS 


The steady and nonsteady pressure coefficients can 
be obtained from the general Bernoulli equation, 


V? dp 
= +f e = F(t) 
2 pp =P 


© = gy + gi + ge is substituted into the Bernoulli 
equation, and terms are dropped in keeping with the 
Keeping only terms 
m/2 results in the 


®, + 


process used to obtain Eq. (1). 
odd with respect to the plane @ = 
approximate “‘lifting’”’ pressure coefficient, 


—2 [(¢e, Uo) Sa (Yo, U,?) | (20) 


C, = 


If normal force coefficient and moment coefficient are 


defined as 


m & "L T 
C*.=— | / C,R cos 6 dé dx 
S 0 0 
-% ft fe 
Cn => | f C, R(xo — x) cos 0 dé dx 
SL 0 0 


because C, is proportional to cos 6, it follows that 


es ze. 
( ) R dx 
S 0 cos @ 


and 


C,=- 


If Eq. (23) is substituted into Eq. (14), there results 


AERONAUTICAL 


SCIENCES—AUGUST, 1951 


Xo T "L c. ) 
Ow = Ci Rx 1x 
Lt SL I Se ” ais 


we are now in a position to obtain the lift and pitching 
moment coefficients for the four motions discussed in 
Section (IV). 


and 


99 


THE Low-FREQUENCY OSCILLATIONS AND STEADY 
MOoTICN 


(VI) 


The task of obtaining solutions for the nonsteady 
supersonic flow about bodies of revolution using Eq. 
(14) is simplified for the case of low-frequency oscilla- 
tions. It can be shown by expanding Eq. (14) into 
powers of reduced frequency that, for harmonic fre- 
quencies of the order of c8/L, the nonsteady potential 
can be represented by the first-order terms in & with an 
error of about 5 per cent or less. For bodies of the 
order of 20 to 30 ft. in length in the Mach Number range 
of 1.5 to 2.5, the first-order terms in k in Eq. (14) would 
represent the nonsteady potential for harmonic fre- 
quencies of the order of 35 to 115 rad. per sec. within 
about 5 per cent. Before proceeding to this expansion 
process, circumvent the improper integrals in Eq. (14 
by introducing the dummy variable u through the 


substitution 


cosh u = (x — &)/Br ( 


va _ . { —iwMr | wr. . 
¢g2 = —cos A exp jot; | B f'(« — Br cosh 1) exp - cosh u f cos sinh «} cosh u du + 
cosh ~'(x/Br) ( cB Cl 
o {° ; ; {—iwMr . . for . 
f(x — Br cosh uw) exp - cosh u ¢ sin sinh w } sinh « du + 
cB cosh ~1(x/fr) ( cB cB 


*0 
¢2 = —cos 6 exp } il | E / j 
cosh ~'(x/Br) 


tw M ro : {—twMr | wr , 
f(x — Br cosh 1) exp - cosh u ¢ cos sinh uw) cosh udu} (24 
cp cosh ~!(x/pr) ] cB cB 


Carrying out the expansion into powers of k, we obtain 


f’(x — Br cosh u) cosh u du — 1MBk X 


*0 
I. \(x/Br) 


If k = Oin Eq. (25), the well-known steady potential results, 


(a) Harmonic Pitching About the Point x) 


f’(x — Br cosh 1) cosh u du 


*0 
—£ cos 6 / 3 
cosh ~!(x/Br) 


Eqs. (25) and (26) will be used to obtain the four nonsteady and steady solutions with which we are concerned. 


(26) 


The equation for the distribution function f(£) is obtained from Eqs. (18a) and (19), 


f(&) = [.S(£) T | [ Uoao + iw (x — 


Eq. (27) is substituted into Eq. (25) and Eq. (25) is substituted into Eq. (20) to obtain the pressure coefficient on 


the body 


Xo) ao | 


f'(x — Br cosh u) du + o¢e | (25) 


(27) 


an¢ 


res] 


and 


The 


and 


Eq 


' Sul 
ing e 
and (: 


bo 
bo 


hing 
“din 


EADY 


eady 
Eq. 
“illa- 
into 
fre- 
ntial 
h an 
the 
ange 
ould 
fre- 
thin 
sion 
(14 
the 


(25) 


(26) 


1ed. 


(27) 


t on 





NONSTEADY SUPERSONIC FLOW 509 


= 2 cos 6 eae) ol; aB Ols aB al2R Ol; 
Cc _ = 3B — t: hs 9 _ & . o 
¥ (« l 0 Ox ° l fh Ox . l f fh l f Ox Tr seal | (<9) 


The integrals /,, appearing in Eq. (28) are presented in the Appendix. These integrals will also appear in the re- 
maining expressions for the steady and nonsteady pressure coefficients to be determined. Eq. (28) is now sub- 
stituted into Eqs. (21) and (22) to obtain the desired derivatives, 


—2 ft lL Ol» Ol; ol, 
Crna = < 28 (1 + ) — M?R — BXxo R dx (29 
SL x '" 20x Ox Ox 
Xo 2 “L 1 Ol Ol; ol; 
Coa = Cas + > ia / | 23 (1 + ‘) — \?R = BXo Rx lx (30 
L SL? Jo tT 20k Ox ow 0) 
The slender body theory expression for pressure co- and 
efficient is obtained by substituting Eq. (18a) into Eq. 
(19) and integrating with respect to R to obtain C,,, = > C.. + 
(v)p—+o- The resulting expression is then combined L P 
with Eqs. (20), (21), and (22) to obtain = {° (1 i. a) eR dx (40) 
‘ bo > 0 Ox Ox 
; S(L) Xo vol. 
Cra = 2 l~ + — (31) 2 ; 
S x LS(L) | rhe corresponding slender body theory expressions 
are 
and 
; ; ’ . ne = 2ZLSL)/S — (Xo/L 
Cua = —2[S(L)/S] [1 — (xo/L)]? (32) Ong [S(L)/S] [1 — (v0/L)] (41) 
: F . ; ’ and 
b) Harmonic Normal Oscillations in Pitch 
From Eqs. (18b) and (19), we have S(L) xo\2 
; Ca, = —2— (1 - ) + 
f(é) = [S(E)/x] wo (33) S L 
Eq. (33) is used in Eqs. (25), (20), (21), and (22), “e 
; S(x) x dx 
respectively, to yield Xo vol. 0 12) 
- (42 
—~9 FL ol; L Sih) L2SXL) 
Cw. == / (a1 — MR ) Ras (34) 
SL 0 On d) Steady Angle of Attack 
and From Eqs. (18d) and (19), 
Xo 2 “L ol / & Y/& y 9 
Cua *§ Cag = (a1 — \PR ) Rx dx (35) f(§) = [S(§)/w] Vow (43) 
lL. SL* Je Ov 
a ' Then, from Eqs. (43), (26), (20), (21), and (22), 
rhe corresponding slender body theory expressions are 
_9 *L 
Cruz = 2(vol./SL) (36) Cra = = (*) R dx (44) 
S o \Ox 
and 
. Xo 28 “L ol 
. S(x) xd Cna = I Cus + LS | (= ') xR dx (45) 
oe) wade ; : 
C . S(L) Yo vol. ws @ (37 ‘4 . 0 - 
mp — < * ov 7° o6) 
5 L S(L)L L*S(L) The slender body theory results are 
c) Steady Pitching About the Point x Cre = 2(S(L)/S] (46) 
Eqs. (18e) and (19) yield and 
f(€) = [S(é)/x] a(é — xo) (38) S(L) Tx vol. J 
Ce = 2 . . 7s l (47) 
§ i £8Ch) 


Substitution of Eq. (38) into Eq. (26) and the result- 

ing equation into Eq. (20) and thence to Eqs. (21) aha j : 

This completes the list of steady and nonsteady 

stability derivatives. These equations will be demon- 

—26 ft (1 Ol» or) Rdx (39) strated in the next section for two different bodies of 
( 


2 - Xo 7: 
hs revolution. 


and (22) gives 


ox Ox 


) 
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(2k) ms iid 
Wich ot 
2,9=0° 
\ 07 b 
"| 
\ ane 
y w(x,t) 
M>I x 
. | wit) “af 
FIG. CYLINDRICAL COORDINATES & ORIENTATION 


of BODY of REVOLUTION 








° 


A,6=0 


hee a a 


wix,f) “Usted + juolx-xeltee 
FIG. 2a PERIODIC PITCHING ABOUT xX, 


X 


= en 


wit) = woe" 
FIG. 2b PERIODIC NORMAL OSCILLATIONS 


x 














A = CONST. 





W(X) = oleUot Q(x ~ Xe) ¢ 
FIG. 2c STEADY PITCHING ABOUT xX, 


ac, x 
e Bo res 
wo\\y = ofUs lH Use 


FIG. 2d STEADY ANGLE OF ATTACK 








RESULTS FOR Two BopIEs OF 


REVOLUTION 


(VII) EXAMPLE 


The equations obtained in the previous section will 
The 


stability derivatives for body (a) of Fig. | will be ob 


now be applied to the two bodies shown in Fig. 1. 


tained using both the first order in frequency equations 
and the slender body theory. Following this, it will 
be shown how modifying the shape of body (a) to body 
(b) changes the stability coefficients. 

Body (a) of Fig. | is generated by revolving a para 
bolic are about the » axis. The equation of this para 
bolic are is 


R = 0.02 (10x 7) 
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Body (b) incorporates the nose of body (a) with a short 


cylindrical section and a slightly ‘‘flaring’’ rear see. 


tion. The equations for the locus of the surface of body 
(b) are 

R = 0.02(10x — x?) (OS x 5) (49a) 

R = 0.50 (5 Sx S 6) (49b) 

R = 0.08715x — 0.0229 (6 Sx FS) (4% 

The point of rotation xo for the steady and nonsteady 

pitching derivative was taken to be x» = 4.78 for both 

bodies. This point was chosen as being typical for 


certain rocket-propelled missiles. 

Fig. 3 presents the values of C,,, and C,,,. for body (a 
using Eqs. (29), (30), (31), and (32). The limitations 
of the slender body theory in predicting the magnitude 
and variation with Mach Number of these stability 
coefficients is apparent. This depends, of course, on 
the validity of the first order in frequency theory. 

Figs. 4, mis GC 
Ca Cry and C,,, for body (a) using both theories. 


5, and 6 present values of C,,;, ( os 
The remarks above relative to the two theories also 
apply to these results. 

An interesting result of these calculations is obtained 
by comparing Figs. 3 and 6, which present the non- 
steady and steady pitching derivatives, respectively. 
In dynamic stability analyses the steady pitching deriv- 
and C,,, are often used rather than the 
If this procedure 


atives C,,, 
nonsteady derivatives C,,, and C,,,. 
were followed for body (a), both the magnitudes of the 
derivatives and the variation of the pitching moment 
derivative with Mach Number would be in error in the 
Mach Number range of Figs. 3 In both cases, 
however, the sign of the pitching moment derivative 
is such as to indicate a damping moment due to a pitch- 


and 6. 


ing angular velocity. 

The slender body theory derivatives for both bodies 
are given in Table | and indicate qualitatively the ef- 
fects of changing the shape of body (a) to that of body 
the same reference 


(b). <All coefficients are based on 


area—the cross-sectional area at x = 5. 


{ 


| 


Table | shows that body (b) has increased lift over 


that of body (a). It is also apparent from the values 
of Cy, that the static center of pressure has moved rear- 


ward with changing body (a) to body (b). Thus, body 
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(b) has an improved margin of static stability over that 
of body (a). 


Both bodies possess steady and nonsteady damping 
in pitching about x». However, body (b) shows an in- 
crease in both C,,,, and C,,,, over those values appropriate 
to body (a). Thus, body (b) has greater damping in 


pitch than body (a). 


The above discussion indicates briefly how the aero- 
dynamic properties of bodies of revolution can be modi- 
fied by varying the shape of the body. In general, it 
can be stated that increasing the base area of a body of 
revolution will improve its steady and nonsteady lift 
and pitching moment characteristics in the supersonic 
range of Mach Numbers. 
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Appendix 

THE INTEGRALS 7, (n = 1, 2, 3, 4) APPEARING IN THE 


STEADY AND NONSTEADY PRESSURE COEFFICIENTS 


Certain integrals appear in any or all of the expres 


sions for pressure coefficient derived in this paper. 


These integrals are: 


TABLE 1 
Slender Body Theory Stability Derivatives 


Body cS Ca G ts 
- a a q 

(a) 0.82 0.93 0.33 —() 

(b) 3.63 0.15 1.47 —0 


(Continued on pu ge 4?) 
Ona Cm Cniy Cm jj 
1.58 —0.13 1.25 —0.002 
3.08 —0.59 1.61 —0.114 








Determination of Rocket-Motor 
Heat-Transfer Coefficients 


by the Trans 


ient Method’ 


STANLEY GREENFIELD? 


North American Aviation, Ine. 


SUMMARY 


An experimental procedure has been developed to measure gas- 
film heat-transfer coefficients in rocket-motor combustion cham- 
bers and nozzles. The procedure utilizes the transient tem- 
perature rise of the uncooled segmented walls of a rocket motor 
to determine the rate of heat transfer. 

This ‘‘transient method”’ was used to determine axial variation 
of film heat-transter coefficients of gases at temperatures up to 
3,000°F. flowing through a de Laval-type supersonic nozzle at 
rates typical of rocket-motor operation. Experimentation re- 
sults of 107 tests are presented 

Forty-nine tests were conducted with high-pressure air at 
1,400°F. and at mass velocities varying between 1.4 X 10° and 
6.8 X 10° lbs. per ft.2 hour. The remaining 58 tests were con- 
ducted on an oxygen-ethyl alcohol rocket motor of 1,000-Ib. 
nominal thrust by passing the combustion gases through the 
test nozzle at flow rates of between 3.5 X 10° and 11.3 X* 105 
Ibs. per ft.2 hour and at estimated temperatures of 3,000°F. 

Experimental values of nozzle gas-film heat-transfer coefficients 
both for the rocket-motor and the lower temperature tests were 
satisfied by the relation 


9 


Ng = 0.029 (G%8/D®-?) c,u%? 


NOMENCLATURE 


A = area (inside surface of nozzle segment), sq.ft 

D = diameter (inside nozzle at segment axial midpoint), ft. 

G = mass velocity (at midpoint of segment), lbs. per ft.* hour 

¢, = specific heat, B.t.u. per lb., °F. 

h = film heat-transfer coefficient, B.t.u. per ft.2 hour, °F. 

qg = rate of heat flux per unit area, B.t.u. per ft.? hour 

t = temperature, °F. 

0 = time; hours 

u = Viscosity, combustion gases, lbs. per ft. hour 

L* = characteristic chamber length; volume of the com- 
bustion chamber to the nozzle throat divided by the 
cross-section area of the throat, in. 

Subscripts 

g = gas 

w = wall, inner segment 

a = average, segment 


INTRODUCTION 


AS THE FIELD OF ROCKET-MOTOR DESIGN advances 
from rule-of-thumb to a more scientific quanti- 


Presented at the 1950 Heat Transfer and Fluid Mechanics 
Institute, Los Angeles, June 28-30, 1950. 

*The author wishes to express his thanks to Joseph Friedman, 
formerly of North American Aviation, Inc., for his assistance and 
guidance in the execution of this papet 


Tt Propulsion Specialist, Aerophysics Laboratory 


tative basis of calculation, the need is apparent for 
reliable engineering values of axial variation of gas- 
film heat-transfer coefficients over the motor internal 
surfaces. 

Although several methods have been proposed for 
such measurements, two procedures appear to be used 
most frequently. The first of these requires careful 
measurement of the steady-state temperature rise of 
the motor coolant as it flows around a given area of 
motor wall. To observe axial variation of heat trans- 
fer, the chambers and nozzles are constructed in several 
short axial sections that are individually cooled. How- 
ever, the increasingly difficult task of measuring pre- 
cisely the small coolant temperature rise and our pre 
liminary knowledge of the coolant-film heat-transfer 
coefficient limit this method as a means for calculating 
combustion gas-film coefficients. 


The second procedure uses the difference of tem- 
perature at two or more precise points within the wall 
of the motor at the same instant of time, from which 
the heat-transfer rate and the gas-film heat-transfer 
coefficient be calculated. The shortcomings of 
measuring accurately the temperature at a given point 
in a metal wall are well known in the field of pyrom- 


may 


etry, and achievement of this goal represents a cardinal 
task. 

A method of measuring the gas-film heat-transfer 
coefficient suggested by the work of London, Nottage, 
and Boelter® appeared to possess substantial advantage 
in its simplicity and direct approach. Briefly, these 
authors utilized the transient temperature decrease of a 
thermal capacitor discharging into a fluid stream to 
determine gas-film coefficients with results that were in 


. . . j 
substantial agreement with those obtained by means 


of other orthodox steady-state procedure. 


Following this principle of the transient method, @ 
procedure has been developed and a copper rocket- 
motor nozzle has been designed and fabricated for gas- 
film coefficient studies in a 1,000-Ib. thrust rocket 
motor. The construction of the nozzle was such that 
it permitted measurements in five axial sections, and, 
in addition, with minor alterations, two to four circum: 
ferential positions could be evaluated in each axial! 


section; a maximum of 20 measurements. 
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ROCKET-MOTOR 


Prior to testing the heat-transfer nozzle (H.T.N.) on 
a rocket motor, it was felt desirable to check the vari- 
ous thermocouple installations and attendant instru- 
mentation and to judge the efficacy of the transient 
method under similar yet more controllable flow condi- 
tions than those provided by the rocket motor. There- 
fore, a series of tests was made first on equipment 
that supplied heated, high-velocity air. The results 
of this lower temperature work served as a guide for 
the rocket-motor tests that followed. 


TEST EQUIPMENT 


Heat-Transfer Nozzle (H.T.N. 

The copper nozzle, shown in Fig. 1, was designed and 
fabricated for installation on a liquid oxygen-ethyl 
alcohol, 1,000-Ib. thrust rocket motor to measure gas- 
film heat-transfer coefficients at five positions along 
the axis of the nozzle. 

From Fig. 2, it will be seen that the internal nozzle 
geometry for a 1,000-Ib. thrust motor at 300 Ibs. per sq. 
in. absolute chamber pressure was machined out of a 
bar of electrolytic copper, and five '/ js-in. wide circum- 
ferential grooves were cut from the outside surface to 
within '/, in. of the These 
acted as insulators between five rings of copper, which, 


inside surface. grooves 
in turn, served as thermal capacitors for the heat re- 
jected to the inside surface by the rocket combustion 
gases. 

A study of circumferential variation of the gas-film 
coefficients is possible by cutting these rings into cir- 
cumferential however, the of this 
work covers only the five axial segments described 


segments ; scope 
above. 

The segmented nozzle was installed for test on two 
pieces of equipment which delivered gases at 1,400 
and 5,000°F., respectively. 


Low-Temperature Test Equipment 


Heated, high-pressure air was furnished by the ex- 
haust of a high-pressure kerosene-air burner, which 
ordinarily supplies the working fluid for a jet pump. 
This installation of the H.T.N. is shown in Figs. 3 and 
1, Air, flowing at velocities varying between 1 and 3 
Ibs. per sec. at 80 to 100 Ibs. per sq.in., was burned 
with ANF 58 (jet-engine fuel) to yield gases at a maxi- 
mum temperature of 1,700°F. The flow of 
measured with a calibrated flat-plate orifice, and the 
flow of fuel was measured with a rotatmeter; the gas 
unshielded chromel- 


air was 


temperature was sensed by an 
alumel thermocouple placed at the centerline of the 
first straightening section (see Fig. 3). 

During the first 37 tests, the straightening sections 
(15.7 in. in length and 3.3 in. inside diameter) were 
During the remaining twelve tests, the 
They were fab- 


water-cooled. 
Straightening sections were uncooled. 
ricated from Type 347 stainless steel with !/4-in. thick 


HEAT-TRANSFER COEFFICIENTS 
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Fic. 2. Detailed cross section of the heat-transfer nozzle 


cylindrical walls. No shielding or insulation was in- 


stalled on these uncooled chambers. 


High-Temperature Test Equipment 


A fixed injector and chamber configuration was used 
throughout the 58 rocket-motor tests. The propellants 
used were liquid oxygen and 92.5 per cent ethyl alcohol. 

The injector was chosen for use in these tests from 
previous experience, which showed that average heat 
transfer to the nozzle, reproducible performance be 
tween runs at the same mixture ratio, and smooth 
burning through a wide range of mixture ratios and 


motor L* could be expected. The specific thrust at 
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Fic. 3. Diagrammatic sketch of segmented nozzle installation 
for operation at low temperature 
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Fic. 4. View of equipment for low-temperature tests. 
300 Ib. per sq.in. absolute chamber pressure with this 
injector was usually 85 per cent of the theoretical at 
optimum mixture ratio. 

A conventional rocket-motor chamber, with L* of 
{S in., was used in all of these heat-transfer expert- 
The cylindrical chamber was_ water-cooled 
A nominal chamber 


ments. 
through a helical cooling passage. 
pressure of 300 Ibs. per sq.in. absolute was used in all 
rocket tests. 

The assembly of the rocket motor is shown in Fig. 5 


during Run N15. 
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APPLICATION OF THE TRANSIENT METHOD 


The method of evaluation of heat-transfer coefficients 
as described by London, et al.,* was used in principle 
for the analysis of the data of this report. 

The circular segments of the H.T.N. acted as ther- 
mal capacitors insulated axially from each other by an 
air gap and essentially insulated on the outside sur- 
face by the surrounding air. Heat was transferred 
from the combustion gases through the gas film to the 
inner wall of the segment and then conducted in a 
radial direction through the segment body. The axial 
variation of the heat transferred in an individual seg- 
ment was not considered. The average temperature 
rise of the segment with time reflected the amount of 
heat which was absorbed by the capacitor. All heat 
transfer is assumed by forced convection only. 

The rate equation for flow of the unit heat flux, g, 
through the gas film is 

q=h,(t, — t,) (1 
and the energy balance on the thermal capacitor may 
be expressed as 

g = (W., A) (di,/d6) (2 
Combining the energy balance and rate equation re- 


sults in 


i= 


3 [W.,/A (t, — t,)] (dt, d6) (; 


g 


Description of Terms 


The terms W and A are geometric or physical prop- 
erties of the copper segments which are readily meas- 
ured. These segment constants are listed in Table |. 
The value of the specific heat for copper at the segment 
average temperature was used in calculations of all 
tests. 

The rate of change of segment average temperature 
with respect to time, d/,/d6, is the slope of the tangent 
to the average temperature-time curve as shown in Fig. 
10. These temperatures were sensed with iron-con 
stantan thermocouples and were recorded on appro- 
priate automatic recording instruments. 

In order to minimize the effect of radiation losses 
from the H.T.N. outside surfaces, tangents were taken 
at the earliest time in the test after steady operating 
conditions had been obtained. Calculations showed 
that the effect of radiation loss on the values of the 
heat-transfer coefficients was less than | per cent when 
the tangent data were taken at the 5th sec. of oper- 


TABLE | 
Physical Constants of the Heat-Transfer Nozzle Segments 


Segment Weight, Segment Inner Area 


Segment Number W, Lbs A, Sq.Ft 
I 3.5 0.049 
II +. ] 0.03 
III 3.9 0.034 
IV > 6 0.042 
\ 1.8 0.081 
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ROCKET-MOTOR HEAT 
ation in the rocket tests and at the 20th sec. of oper- 
ation in the low-temperature equipment. 

The thermocouple within each axial segment was 
installed at the location of the segment average tem- 
perature as calculated from Schmidt plots of the noz- 
ale wall. Procedures outlined by Schmidt* and Perry 
and Berggren‘ were followed in the preparation of these 
temperature gradient plots, wherein a scale distortion 
of the abscissa is required to apply the Schmidt method 
to plane radial heat flow. The segments of the H.T.N. 
were considered to be finite hollow cylinders, and Perry 
and Berggren have shown that the use of a logarithmic 
scale in plotting equal finite increments of cylinder 
radius makes possible the graphical application of the 
Schmidt method to cylindrical systems in the same 
fashion as to slabs and with equal precision. 

The area beneath the curve of temperature versus 
radial thickness of the nozzle wall segment was graphi 
cally integrated and averaged for the position of the 
average temperature after the radial thickness coordi 
nate was weighted to account for the variation of the 
mass of metal with respect to thickness. The graphi 
cal integration described, therefore, is essentially a 
summation and an average of the heat content of the 
segment. The heat 
content and temperature at a particular instant of time 


positions of segment average 
are alike. 

For these Schmidt plots the values of specific heat 
and thermal conductivity of copper used were assumed 
constant at values corresponding to the average tem- 
perature experienced by the metal segment during the 
test. The gas-film coefficient was arbitrarily calcu- 
lated from Eq. (4)) of MeAdams,!' and a gas temperature 
of 4,500°F. In the event experimental 
values of the gas-film coefficient radically disagreed 
with this selection, a reiteration process would be ap 
plied on the Schmidt Plot using the experimental 


was used. 


values of /, and ¢,. 

The temperature of the gas, /,, which affects heat 
transfer across the gas film, was a variable factor de- 
pendent upon the test equipment and the position in 
the nozzle. At the three segments in the divergent 
portion of the nozzle, the gas temperatures were con- 
sidered to be equal to that of an adiabatic wall along- 
side the same fluid as the segment wall. The tem- 
perature of the adiabatic wall was computed from the 
main-stream conditions and a recovery factor from a 

given by Kaye, Keenan, McAdams.’ 
heat-transfer coefficient was insensitive to vari- 
the the factor 
because the difference segment 
perature and adiabatic wall temperature was relatively 
It was appreciated, nevertheless, that the 
measured by Kaye, et al., would 
not apply the flow 
supersonic nozzle, however, application of a recovery 


curve and 
The 
ations in 


selected value of recovery 


between wall tem- 
large. 
recovery factors 
directly to conditions in a 
factor was considered necessary and no more correct 


information was available. 
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test 


| a typical rocket-motor heat-transfet 


For the two segments in the convergent portion of 
the nozzle and the throat, the stagnation or total tem 
perature of the gas was taken as/,. At these segments 
the gas velocity was equal to a maximum of only the 
speed of sound and the recovery factor was therefore 


of little significance. 


The value of the main-stream gas temperature in 
the rocket motor, equal to the stagnation temperature 
in the nozzle, can be calculated from thermodynamic 
data when complete chemical equilibrium and uniform 
mixing of the propellant gas components aro assumed. 
Test runs of rocket motors showed, however, that these 
assumptions were not strictly correct and that the 
effective combustion temperature, as judged by rocket 
motor specific thrust, was lower than the calculated 
thermodynamic value. A well-known thermodynamic 
relation showed rocket-motor specific thrust propor 
tional to the square root of the absolute combustion 
temperature. that 
ratio of specific heats and molecular weight were neg 


It was assumed variations of the 
ligible and that the square-root relation of measured 
specific thrust of the rocket motor was used to calculate 
the effective combustion temperature. 


Combustion temperature of the gases from the low 
temperature equipment was sensed with an unshielded 
chromel-alumel thermocouple. Corrections were ap 
plied to the measured gas temperature to account for 
loss and factor at the thermo 


radiation recovery 


couple. 


For the rocket-motor tests, temperatures of the inner 
wall of the H.T.N., ¢,, adjacent to the combustion 
gases at each axial segment were estimated from the 
Schmidt plots. 
inner wall and the position of the segment average was 


The temperature difference between the 


obtained by multiplying the ratio of experimental gas 
temperature used in the Schmidt plots, 4,500°F., by 
the inner wall-to-segment average temperature drop as 
given by the Schmidt plot at 5 sec. time. This ap 
proximation was considered to be valid within the 
accuracy of the results, since the experimental and 
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assumed values of the gas-film coefficient checked 
closely. 

Wall temperatures in the nozzle during the low- 
temperature ram-jet equipment tests were extrapolated 
linearly from the temperature drop that was recorded 
between two thermocouples at different depths in the 
throat segment. One temperature was taken at a 
depth of '/2 in. to measure average segment temper- 
ature, and another was taken 1 in. in from the outside 
surface. The validity of this linear extrapolation of 
the temperature at the wall was supported by the low 
values of the gas-film coefficient and the gas tempera- 
ture during these tests. 


DISCUSSION OF RESULTS 


Low-Temperature Tests 

It will be recalled these tests were conducted to 
achieve a degree of familiarity with the heat-transfer 
equipment and to evaluate better the effectiveness of 
the transient method. The low-temperature equip- 
ment provided control in variation of temperature 
and flow, conditions not easily obtained in a rocket 


motor. 


The complete series was conducted at a gas tem- 
perature of 1,400°F. (total temperature) and a gas 
flow rate varying between 1 and 35 Ibs. per sec. Two 
factors were considered of importance in selecting these 
conditions: (1) The equipment operators could reach 
steady burning conditions most quickly, and (2) the 
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gas Reynolds Numbers in the nozzle were of the same 
order of magnitude as those in rocket tests. 


Results are shown plotted in Figs. 6 and 7 in the form 
of Nusselt and Reynolds Numbers based on inside 
diameter at the axial midpoint of each segment as the 
characteristic dimension. The physical properties of 
the gas were assumed to be essentially that of air at 
1,400°F., since the amount of fuel burned with air to 
form the hot combustion gases was far below stoichi- 
ometric and represented less than 5 per cent of the 
total gas flow. 
corresponding to cooled and uncooled upstream chamber 
conditions. No correction is applied to account for 
the maximum 3 per cent heat loss from the gases to the 


walls of the straightening sections. 


The inlet chamber conditions for the data of Fig. 6 
closely approached that of a cooled rocket motor, and 
it was noted that a straight-line correlation of 85 per 
cent of the experimental points was obtained within 
+15 per cent. This scatter of the experimental data 
could be anticipated by considering the nonuniform 
flow field at the inlet to the heat-transfer nozzle. After 
passing through a 90° elbow and a 6- to 3.3-in. inside 
diameter transition section, the 4.7 diameter lengths 
of straightening section could not be expected to estab- 
lish , reproducible gas velocity and temperature pro- 
files at the nozzle entrance. Additional calming lengths 
would have been desired but could not be installed be- 


cause of space limitation. 
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The data is presented in two groups | 
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The conditions in an uncooled rocket motor were 
simulated in a second series of twelve tests when the 
straightening sections were uncooled. The data of 
these runs are plotted as Nusselt Numbers vs. Reynolds 
Numbers in Fig. 7 showing line B-B as in Fig. 6. The 
values of Nusselt Number show a trend parallel to, but 
about 50 per cent higher than, the line B-B. Inspec- 
tion of the Nusselt Numbers of the individual nozzle 
segments showed the highest values for Segment I 
(see Fig. 2), with progressively lower values at Seg- 
ments II and III and, finally, values corresponding to 
B-B for Segments IV and V. The higher trend in 
Nusselt Numbers is attributed to the nature of the 
nozzle inlet gas temperature profile during this test 


series. The stainless-steel metal surfaces of the un- 


cooled straightening section adjacent to the gases 
quickly reached an elevated temperature causing the 
combustion gases in and near the wall boundary layer to 
reach a higher average temperature as compared to 
previous tests with cooled straightening sections. In 
effect, a discontinuity of metal surface temperature 
occurred at the entrance of the nozzle where lower 
surface temperatures were expected on the highly 
conductive copper. The result would be a higher ap- 
parent value of Nusselt Number at the leading seg- 
ment with subsequent decrease in the effect as the gas 
temperature profile became adjusted to the new bound- 
ary conditions in the nozzle. 

The correlation of data, particularly that taken with 
cooled straightening sections, and the reproducibility 
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of test results reflected satisfactory operation of the 
segmented copper nozzle with its attendant instru- 
mentation and furthermore asserted the potentiality 
of the transient method for measurement of gas-film 
heat-transfer coefficients in a rocket motor. 


High-Temperature Rocket-Motor Tests 


Correlation of the rocket-motor heat-transfer data 
was attempted through plots of dimensionless group- 
ings, Nusselt and Reynolds Numbers, based on inside 
diameter of the nozzle, and length along the nozzle 
wall, as the characteristic dimension, D. The points 
fell in a nonuniform fashion, however, and it was not 
possible to fit a representative continuous curve or line 
to the result. Probable source of this poor correlation 
lies in the lack of accurate values for the thermal con- 
ductivity and viscosity of the combustion gases. At 
the time of this report, the author is unaware of re 
liable extrapolation of transport properties of gas mix- 
tures or experimental measurements of conductivity 
and viscosity at temperatures as high as 3,000°F. 

Because of their inaccuracy, the use of dimensionless 
groups was set aside for the present in favor of a more 
elementary plot. A straight-line relation of data was 
obtained by plotting the segment gas-film coefficients 
versus the ratio G°-5/D°-? as shown in Fig. 8, with D 
the inside diameter at the midpoint of the segment. 
The validity of such a plot was further demonstrated 
by the correlation of the rocket data with the data from 
the low-temperature tests at 1,400°F. (see Fig. 9). 


AUGUST, 1951 


The specific heat and the viscosity of the combustion 
gases at the recovery temperature were applied to the 
ratio G®-* D°-? of Fig. 8 in order to adjust for the dif- 
ference in temperature of the two sets of data. The 
gas viscosity entered only to the two-tenths power 
for this application and the effect of inaccurate values 
was minimized. The specific heat of the gas mixture 
at the higher temperature was readily calculable and 
probably accurate to within 5 per cent. 
The equation of the straight line in Fig. 9 is 


h, = 0.029 (GP, D®) ete? ! 


and it is noted that 75 per cent of the rocket-motor 
test points are satisfied by the above equation with 
+15 percent. The data generally fell into five groups, 
numbered in Fig. 9, corresponding to the five segments 
of the nozzle (see Fig. 2). 

Heat losses from an individual segment radially to 
the atmosphere and axially to an adjoining segment 
through the -in. web at the inside wall have been 
shown through preliminary calculations to be less 
than 2 per cent and were neglected in this report. 
There is an additional heat loss through forced con 
vection caused by a movement of air over the nozzle 
induced by the action of the jet flame. Experimental 
checks of the magnitude of this convection loss have 
not been made at this time, since initial estimates 
showed an effect less than 2 percent. The total of the 
error in the values of gas-film heat-transfer coefficients 
as a result of the above three heat losses was estimated 
to be not more than 5 per cent. However, the con 
sistently low values of the gas-film coefficient in the 
exit Segment V have been in part attributed to the ex 
cessive amount of turbulence over the outer and end 
surfaces of this particular portion of the nozzle. There 
fore, an experimental evaluation of the convection loss 
is planned at an early date to corroborate the prelim 
inary calculations. 

The less than 5 per cent decrease of gas temperature 
due to transfer of heat to the cooled rocket chamber 
walls is neglected. The effect of change in rocket 
motor mixture ratio on the nozzle gas-film coefficients 
has not been emphasized, since, at the present time, 
information for all the mixture ratios investigated is 
not complete enough to make definite conclusions. 

The information presented here illustrates an initial 
phase of a series of rocket-motor heat-transfer measure- 
ments that are anticipated with the oxygen-alcohol 
propellant system using the transient method. 

It is appreciated that application of the data reported 
here to rocket-motor-nozzle heat transfer in general is 
limited, since a fixed injector and L* configuration were 
used throughout. 

When further test data are recorded with a variety 
of injectors and motor L*, a more complete picture of 
the heat-transfer field will be available to rocket motor 
designers. 


(Continued on page 526) 
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Aerodynamics of a Supersonic Rectangular 
Wing Striking a Sharp-Edged Gust 


THEODORE R. GOODMAN* 


Cornell Aeronautical Laboratory, Ine. 


ABSTRACT 
The method of C. S. Gardner for solving nonsteady finite wing 
problems at supersonic speeds is briefly described. An expression 
for the pressure is derived in terms of the velocity potential in 
the transformed variables used by Gardner. Nonsteady wing 
problems are solved, according to Gardner, by solving two three 
dimensional boundary-value problems. In the case of the sharp 
edged gust, the first problem is a symmetrical finite wing problem, 
and the second is an antisymmetrical problem based on the solu 
tion to the first problem. Expressions are derived for the pres 
sure in the various regions on the wing, and these are shown to 
be conical. Finally, the pressures are integrated, and formulas 
are given for the lift and leading-edge moment as functions of 
time, aspect ratio, and Mach Number 


INTRODUCTION 


eons METHOD PRESENTED RECENTLY by Gardner! 
for solving nonsteady finite wing problems at 
supersonic speeds is particularly adaptable to the wing 
of rectangular plan form. This method, which reduces 
a nonsteady finite wing problem to two “‘steady”’ finite 
wing problems, will be used here to solve for the pres- 
sures on a rectangular wing striking a sharp-edged gust. 

The problem of a sharp-edged gust on an infinite wing 
at supersonic speeds has been solved by Schwarz* and 
others.*-® Moreover, Miles has solved the problem 
fora rectangular wing.’ The method he uses, however, 
permits him to derive only the chordwise lift distribu 
tion. The actual pressure at a point is not given. 
Not only are the pressures required in certain applica 
tions, but physical insight is gained into the nature of the 
behavior of a sharp-edged gust interacting with a wing 
tip when the local pressure at a point is known. 


THE METHOD OF (JARDNER 


Without deriving the method, its essential points will 
be presented here. With rectilinear coordinates x, y, 2 
fixed on the wing (x in the downstream direction and 
y in the spanwise direction), the boundary-value prob- 
lem satisfied by @ (the perturbation velocity potential) 


is, for (x, y) in WW, 


lim @ = lim @, = Uf_(x, y, t) + fx, 9,8) (1) 
om O4 zs 0) 
forx < 0 
Hx, ¥, 3,1) =U (2) 
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(MV? — l der = Pi, — @2z T : Prt : r x2 * 0 (3) 


where IV is the projection of the wing in the xy-plane, 
s = f(x, y, t) describes the surface of the wing, ¢ is the 
velocity of sound, and .V/ is the Mach Number of the 
free stream (JJ = Uc). 

This boundary-value problem involves four inde 
pendent variables. Gardner's method reduces the 
problem to two boundary-value problems, each one of 
which involves only three independent variables. 

First, transform the variables x and ¢ to the variables 
v’ and t’, where 


x’ = x(M? — 1)7' ) 


= (Mx —- (Mt —-pnajare—y-“f 


Let g = Uf, + f; Now drop the primes on the 
transformed variables, vx’ and ¢’. Then Eqs. (1), (2), 
and (3) become, for a rectangular tip, for y > 0, 


lim @. = lim ¢, = g(x, y, ¢) (5) 


z—> 0+ s—> 0 
(where g = 0 for x < 0), and for 
} ¢=0 (6) 
On. — Oy — @: o,, = V (7) 
Gardner's method consists in introducing the variable 


¢ and the quantities y and x. With a slight modifi- 


bs 
cation of Gardner's presentation, x satisfies the follow 
ing boundary-value problem, for — > 0, 


Xrr — Xe Xee 0 (S 

for — = 0, 
x = g(x, y, t) (9 

forx < 0, 
x =0 (10) 


y satisfies the boundary-value problem 


Vie — Yn — vez = 0 (11 
for \ > e 
lin yy, = lim y x(é, x, y, t) (12 
atl beastie 
for —& > x, 
y =0 (13 
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Fic. 1. Wing section entering vertical gust. 
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Fic. 2. Regions of the x - ¢ ‘“‘wing.” 
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Fic. 3. The second problem—part 1. 
The velocity potential is found by setting & = 0, such 
that 
o(x, y, 2, t) = (0, x, y, 2, t) (14) 
The pressure in terms of physical x/ variables is given 
by 
b = p[(0¢/Ot) + U(0d/Ox) | (15) 
In terms of the transformed variables the pressure is 
pu Op Op 
p= ( + M (16) 
MV M? — 1 \ol On 


THE First PROBLEM 
From Fig. 1 it is seen that, in terms of physical x, ¢ 
variables, the boundary condition is 


x < Utl 
x > Uts 


g=v for 


g = 0 for 
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In terms of transformed variables [Eqs. (4) ], 


g=% for x > Mt 
=a : ( (18) 
g=0 for x < Mt 
It is convenient to define a quantity o such that 
Oa /O0E = x (19) 


g satisfies the wave equation [Eq. (8)], and then 


Eq. (9) becomes 
Oa /0E = g(x, y, t) (20 


The boundary-value problem defined by Eqs. (Ss), 
(20), and (10) can be thought of as a “symmetrical” 
supersonic “‘wing’’ problem in x, /, — space, where x 
is in the ‘‘downstream’’ direction, ¢ is in the “span- 
wise”’ direction, and £ is normal to the wing as shown in 
Fig. 2. 

The solution for o can be obtained in terms of simple 


sources 


| ‘ r Vo dx, dt; 
c= = = = (21 
a , (x — 1)? — (t — 4)? -— &] 


where 7 is the shaded region shown in Fig. 2 provided 
x, t, € lies within the ‘‘Mach cone”’ from the apex of the 
“wing.’’ Call this phase (2). If x, t, & lies to the right 
of the Mach cone, then o = 0. If x, ¢, & lies to the 
left of the Mach cone, o must be ‘“‘two-dimensional”’ 











i.e., independent of ¢. Call this phase (1). Ther, 
clearly, 
og. = 0; on Mach cone left (22 
oz = 0 on Mach cone right (23 
where the subscripts refer to the phase. 
g 
cs 
@) men, 
S ea Vx -t? 
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Phase (1) 


Putting in the limits of integration on Eq. (21), there is found, for o, 


Vo i ee i ss vil he + dx, 
aaa dt, : (24) 
TIJt-—-Vxt-# 0 V(x — mu)? —(¢ —4)? -—-& 


Phase (2 
To obtain the proper region 7 in Eq. (19), use the result of integrating Eq. (24) and subtract the extraneous area 


from Fig. 2; then we obtain 


vy {7 7 Mt dx, 
oo = H(— — x) + / dt | + 
T JO 0 V(x — 4)? — (tt -— th)? - & 
Vo "t+ \ di = "x \ f f 7 dx; 
/ dt, | . (26) 
7 T 0 Vv (x — x)" — (¢t — t,)° — gE 


1 = [Mx —t - V(x — M1)? 4+ OP — DEV/OP — 1) (27) 


where 


Performing the x, integration, there is obtained 


oe a x — Mh ty (itVEE ‘ 
o2 = m(§ — xX) — / com-* — dt, + cosh ~ dt, (28) 
' ¥¢G- +? 7 Jo V(it—at+e# 


It will be seen later that, in order to find the pressure on the wing, ¢ itself is not necessary but only the combina- 


tiona, ++ Mo,. We therefore proceed to find this quantity. 


Oo» Oo» Vo er tO re) x — Mt, Vo x — Mr 
+ VM = —My — | + M cosh~! ; dt; — cosh~! 
Ol Ox Ls 9 ot Ox V(t — th)? + # V(t — 1)? + #2 


(° + M =) é Vo i alae (- 41 M >.) h-! . It; (29 
] ] -osh 2¢ 
Ot Ox i T I Ol Ox , Vv (f ron t,)? + £ — 


But 
[(0/ot) + M(0/dx) |r = 1 (30) 


Also 


*t+Vx2-# Q x iV e—F'¢ x 
/ cosh~! dty = — cosh~! dt, (31) 
J of V(iti—-hy?+e2 0 Of; Vit—-ne+e 


The last integral in Eq. (20) is easily evaluated, and so there is obtained 


Oo» Oe vu {" (Mt — x) (t-— hh) + ME] dh 
+M— = - / . 
of On trJo V(x — Mh)? — (t — 4)? — Bf -— 4)? + &] 
Vo x — Mr, Vo x Mv { 
cosh~! + cosh™! _ cos! (32) 
Viti—n?+e 7 Ve+e Vx? — £2 
The substitutions 
t—t, = ~étan@ 


and 
v9 = (x — Mt) cos 6+ WME sin 0 


make the first term readily integrable. Finally, there is obtained 


(Oa2/0t) + M(002/0x) = —(Muo/r) cos! (t/V x? — & (33) 
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From Eq. (25) there is found 
(Oo,/Ot) + A(O0, Ox) = — Mav (34 


9 


On the Mach cone right (¢ = V x? — &), Eq. (23) is satisfied; on the Mach cone left (¢ = —Vx* — &), Eq. 


(22) is satisfied. 


THE SECOND PROBLEM 


There are two parts to the second problem. Part | occurs when / > 0; part 2, when/ < 0. Let us solve part | 
first. 

Part | can be solved with sources by Evvard’s method” when 0 < y < ¥V x? — £2. When y > Vx? — #*, the 
solution becomes two-dimensional—i.e., independent of y. In all cases we shall let s — O so that we find solutions 
on the wing only. We obtain ¢ directly by letting = 0. The solution for y > V x* — f is obtained by integrat 


ing over the entire fore “‘Mach cone” as shown in Fig. 3. 


(o( £1, x, t)dy, 


6 

ll 

| 
| — 
— 2 

— 

rm 
Mee 
4 
Ser || >e 


7 — (y — 9)? 
*a/at—t 
= / x2( £1, x,t) dé, (30 
—¥ © 
Using Eqs. (19) and (23) we obtain 
Qo = — (0, X, t) (.94 
From Eq. (33) we find 
(Om /Ot) + MW(0¢/0x) = (Mvo/7) cos! (t/x) (38 


i] 


Now we obtain the solution for y < V x? — ¢ by using Evvard’s method’ and integrating over the shaded region 
shown in Fig. 3. Again setting = 0, we obtain, for ¢, 
| oe oe 
QoQ = —o»(0, 3) + x2 d 
1 i 


“ti —y dy, 
' (39 
-—e Vr- ¥ - wy 
After some manipulation there is obtained 
Vy {V¥?-® olf, x, t) dé 
@ = —o.(0, x, tf) + — ° (40 
9 # ; 


From Eqs. (33) and (23) we obtain (dropping the subscript on &) 


Str 


Oo Op Mvp t Muy “Vv x8-# eos—! (t/V x? — £*)dé 
+ M = cos! - Vy = (41 
Ol Ox TT x T a4 twte-~y 
We now proceed to part 2 (¢ < 0). 
In part 2, forx > & > Vx? — t,x = x1; whereas in part |, forx > & > Vx? — #, x = 0. The regions in 


duced in the second problem are shown and labeled for convenient reference in Fig. 4. Except for region (1), 


Evvard’s method becomes cumbersome, and for the other regions the writer’s method®: '’ will be used instead. 
Let us first find the solution in region (B). This may be found by simple sources. Setting & = 0, we obtain 


a —| f x2 — 42 dt, f th -x0(&1, x, 1) dy, 7 | f dé f t& x4 (£1, x, L) dy, (49 
r Jo vi VE? —(y—m)? TH Vee 9 Sy-n VE? — (y — 91)? 
After using Eqs. (22) and (25) 
@ = —o(0, x, t) (45 
Therefore 
o, + Mo, = (Mvo/r) cos (t/x) (44 


Therefore, there is no difference between parts (1) and (2) in the two-dimensional region {see Eq. (38) ]. 
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part 2. Diagram for finding ys. 
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Fic. 6. The second problem —part 2. Diagram for finding y, 


The solutions in regions (A) and (1) are those for a 
steady two-dimensional and tip region, respectively, 
since x; is constant. These solutions have been ob- 
tained elsewhere [e.g., Eq. (20) of reference 5] and will 


inerely be restated here. 


Wa = Vo(x — &) 
¥1 [eVvyw«e —e- y) + 
A(x — &) tan“! Vy/(x — = — y)] 


Letting € = Oand differentiating, we find, for ¢, + .\/¢, 


(Od, Ot) + (Od, Ox) = Map (45) 
O¢ O Mvp 2y 
= 4 oe 2 cos (1 = 2) (46) 
Or Ox T x 


We shall now proceed to find the solution in region 
(2) by the method of reference 9. 

Applying the method of reference 9 to the problem 
shown in Fig. 5 we arrive at the double integral equa- 


tion 


> : ; 
JI Jc ro (Wr — v2) diddy = 32 (xr — x2) (47) 


where 


(48) 


>>) 
\| 


— Vi )2 
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After transforming to ‘Mach coordinates’ and soly 
ing each of the double integral equations, we arrive at 


9 ] ~ (0/0) (x1 — x2) dé dy 
~ (hi — 2) = ; 
0) 7 2) V (& — §)*? — (vw — y)? 
(49) 
Here, the roles of &, y; and &, y have changed places. 
Performing the indicated integration and dropping 


the subscript on £, we arrive at 


(Oo O£) (Wy - 2) = 2 = Wi (50) 


or 


Vi — fo = 7 — (51) 


The analysis for finding ¥3 yields the relation ~; = 
yo. Hence, for — = 0, one need not differentiate be- 
9 


tween regions (2) and (3) but can state, for V we P 
a> < HF =, 


Mu t . 
cos”! — ] (52) 
T x 


The analysis for finding Ys leads to the integral equa- 


tion (see Fig. 6) 


/ | K (Ws = Yo) dé, dy, = 
4 
a) K( 2 = 1) dé, dy, (53) 
7) 


After considerable manipulation, the solution to this 


QD; + Mo, = 


equation is found to be 


VY Vx? i (ao, — ox) dé, 
i Y= 3 
w t+) e-fva— f= ¥ 


After letting & = O and using Eq. (33), we find 


(54) 


Od; 1 Ods _ 


Muo | t 
M =— cos 
Ot Ox T x 


+ 


. 272 » ¥9 » 
ve" cos7! (t/V x? — &) dé es 
(ded) 
y EN - — y 


Ss 


An analysis of region (5) shows that ¥; = ¥s. Hence, 
the Mach line bordering regions (2) and (3) and (4) and 
(5) does not exist as a demarcation of solutions. Its 
reflection from the line y = 0 cannot exist either, and 
Fig. 4 should be changed. This changed diagram 1s 


shown in Fig. 7. 
THE PHYSICAL WING 


We shall now recapitulate the multitudinous solutions 
obtained and interpret them in terms of regions on the 


physical wing. The line y = Vx* — f° becomes, in 


terms of physical x, ¢ variables, 
(56) 


(x — Ut)? + y? = (ct)? 
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Fic. 7. The second problem—part 2. Extraneous Mach lines Fic. 8. Regions on physical wing. 


removed. 


This defines an acoustic circle within which the solu- 
tion is three-dimensional and time dependent; here, 
the effect of the wing tip and gust edge is completely felt. 

The regions on the physical wing are depicted in Fig. 
8. The lines are labeled with the physical x, ¢ variables 
and B = V M?- 1. 

As time progresses, the lines (U’ — c)t, Ut, (U + c)t 
move downstream on the wing; the circle moves down- 
stream also, and its radius increases in such a manner 
that the Mach cone from the origin is generated tangent 
to the resulting family of circles. In regions R; and 


R; steady flow has been established. In region R, 


(45) Pi = pUvp B 
(38), (44) po = (pU0/7B) cos (t/x) 
(46) ps = (pU%/mB) cos [1 — (2y/x)] 
? pv 2y t 
(52) fA = cos~! {1 — —] — cos! [{ — 
Lae) £ x 
~ pv 2y f 
(55) Ps = cos~' {1 — — — cos”! : 
re x x 
Ux t y “V s*~F ons! 
(41) ps = : eos i te ee / 
Tk Y Tv y 


Examination of the pressures in the various regions 
Note that the 
Note 


discloses the character of the solution. 
pressure is continuous from region to region. 
also that the pressure reduces to zero along the line 
i.e., it 1S 


y = 0. The pressure is conical in all regions 


a function of x/fand y tonly. In order to gain greater 
insight into the behavior of the pressures, they have 
been determined numerically for JJ = 2, and a diagram- 


matic drawing of the loading is given in Fig. 9. 


the flow is nonsteady but two-dimensional. In regions 
R; and R¢ the flow is three-dimensional nonsteady. In 
region R, the flow consists of two parts—one is three- 
dimensional steady; the other is two-dimensional non- 
steady. These considerations have been arrived at 
by considering the pressures in each region. 
(16) and the various analyses throughout the paper, 
one arrives at the following list of the pressures: Here, 
the subscript number coincides with the numbered re- 
the number preceding each 


From Eq. 


gions as depicted on Fig. 8; 
pressure is the equation number from which the pres- 
sure was taken; x and ¢ are transformed variables. 


(59) 


(60 


y OV x? 2 eos! (\t Vx — ¢?) dé 
= (61 
J» EVE-—y 


THE LirT AND MOMENT 


Although the pressures in regions (5) and (6) cannot 
be expressed in terms of tabulated functions, if the 
pressures are integrated in the spanwise (y) direction, 
the resulting spanwise lift can be expressed in terms of 


elementary functions. This will now be demonstrated. 


The Spanwise Lift 
Refer to Fig. 8 


Let L, represent 


There are three cases to consider. 
and use the transformed x, ¢ variables. 
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AERODYNAMICS OF A SUPERSONIC RECTANGULAR WING 


the spanwise lift, and let b/2 be the semispan. Then, 9, 


case (1) occurs when ¢ < —vwx, and a 


eee. eee 
a b/2 
L;, = 4 (f. ps dy + f pi dy) (63) me 

0 x 








Case (2) occurs when —x < ¢ < Oand \ 
V x?-2 ” Ma 

L= 4 ps dy + prdy + 

i 0 Vx?-P 

*b/2 
J pe dy) (64) 
x \ 
Case (3) occurs when x > ¢ > Oand 
rT 
x? 22 *b/2 : ve | as 
L, = 4 pe dy + p2 dy (65) 4 > 
0 V x?*-le 2 

1 he factor 4 appears in each ol these expressions for gyg_ 7° 3 ae 
the following reasons: A factor of 2 includes both upper q 
and lower surfaces of the wing, and another factor of 
2 accounts for the other semispan of the wing. The . A 

: ° : . . . t 

subscripts on L, refer to the case. Substituting Eqs. . 
aie 9\ 2 > 6 2p . “¢ Fic. 9. Pressure distribution in perspective. 
(57)-(62) into Eqs. (63)—(65) leads to the spanwise lift. 2 : _— 
In order to avoid difficulties with branch points and Hence, it is not necessary to distinguish between 
multivalued functions, it is necessary when integrating L. and L,.. L, is equivalent to the y that Miles’ ob- 


ps to use the absolute value of ¢ in the expression for /s tains and is proportional to it. In order to obtain the 
as given in Eq. (61). The results of the integration are ift and leading-edge moment, first transform back to 
physical x, ¢ variables and then integrate in the chord- 


L,, = (2pU%/ 8) (6 — x) (66) : sine . ; 
wise direction with a lever arm x (physical) for the 
2pUv (b t t x nt : 7; » lever ; ; itt The 
oe ee pury ( cos? - + - — ) (67) moment and without the lever arm for the lift. The 
B 1 x 2 2 results for the lift are 
i Br l Br M 
= 1 — 0O<r< (68) 
Lo M 2 MBA.R.:) M+ 1 


L | - M? — (M2 — 1)r Br [nx a M(1 — 7) 
j.« ( M wer E sis , * 


l (M — 1) (MM? — 1)r? 2(M? — 1)r M ' M 
— + (\J — 1) ( = ee (69) 
{BA.R. MP, M M+1 M-— 1 


L/Ly = 1 — (1/28A.R.) ((M/M—1)<7r< @] (70) 


where Ly (= 2prpbl/B) is the lift on a two-dimensional wing at angle of attack u/U,7 = Ut/l and / is the chord. 


The results for the leading-edge moment are 


Mur Br- 2Br M 
= (1 — ox ¢< (71) 
My M 3MBA.R. M+ 1 
Mit | = — (IP? - ve) Br ofa 70) 24 Br? ; M(r — 1) 
= ‘OS — -~— jf- — é ‘OS 
i «= | \ M _ Ce a oe , 
l 73(M — 1)3(M + 1) 3(M? — 1)r M M 2 
= + 2(M — 1) <r<— (72) 
6BA.R. M: M M+ 1 M-1 
Mir/Mo = 1 — (2/3BA.R.) [M/(M —1)<7< @] (73) 


where \/, (= pl’wb/?/B) is the leading edge moment on a two-dimensional wing at angle of attack v/ U’. 


These expressions for lift and moment agree with those obtained by Miles.’ 
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DISCUSSION 


The solution presented here is valid within the as- 
sumptions of linear theory. The external regions must 
not interact, however, and this condition will be met 


provided BA.R. > 1. 
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Determination of Rocket-Motor Heat-Transfer Coefficients 
by the Transient. Method 


(Continued from page 518) 


CONCLUSIONS 


(1) Using the rocket-motor injector, motor L*, 
and propellant combination reported in this work, 
nozzle gas-film heat-transfer coefficients can be pre- 
dicted within the accuracy required for engineering 
estimates (4). Coincidentally, the 
cepted equation factor for straight pipe heat transfer 


from Eq. ac- 
as given in McAdams! is 0.027 in place of 0.029 used 
above. 

(2) Until 
become available for the viscosity and thermal con- 


more accurate high-temperature data 


ductivity of rocket-motor combustion gas mixtures, 
satisfactory correlation of rocket-motor heat-transfer 
data may be achieved by the use of specific heat data 
in the above type of equation. 

(3) The 
measured under the controlled low-temperature con- 


correlation of nozzle heat-transfer data 


ditions of this report reflected satisfactory operation of 
the segmented copper nozzle with its attendant in- 


strumentation and has demonstrated the potentiality 
of the transient method for measurement of rocket- 


motor gas-film heat-transfer coefficients. 
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Similarity Laws for Slender Bodies of 
Revolution in Hypersonic Flows 


STANFORD E. NEICE* ann DORRIS M. EHRET?T 
Ames Aeronautical Laboratory, NoA.CA. 


ABSTRAC1 


The hypersonic similarity law, as developed by Tsien for 
bodies of revolution at zero angle of attack, is tested by com 
paring pressure distributions obtained using the method of 
characteristics. The results indicate that, in spite of the assump 
tions made in its development, the hypersonic similarity law is 
valid for relatively wide ranges of Mach Number and fineness 
Similarity laws for hypersonic flow about slender pointed 
Both 


the potential flow and the viscous effects of the transverse com 


ratio 
bodies of revolution at angle of attack are also developed. 


ponent of the flow are considered in this development 


INTRODUCTION 


_— LAWS ARE IMPORTANT to the aerody- 
namicist because they provide a useful guide for 
the correlation and extrapolation of data. These 


laws, for example, show how data for a given body, at a 
given Mach Number, can be applied to a related body 
at another Mach Number within the same Mach Num- 
ber range. Similarity laws of this type have been de- 
veloped for the different Mach Number ranges. In 
subsonic flow there is, for example, the familiar Prandtl- 
Glauert rule. A similar rule is available for low super- 
sonic flow. The validity of these laws has been es- 
tablished through repeated checks with theory and ex- 
periment. More recently, similarity laws have been 
developed for transonic and hypersonic flows. The 
validity of the transonic law for two-dimensional flow is 
now being investigated by Liepmann and others. 

The hypersonic similarity law is currently under in- 
vestigation at the Ames Laboratory of the N.A.C.A. 
The purpose of this investigation is to evaluate the ac- 
curacy of the law for nonlifting bodies of revolution and 
to determine the form of the law for more general three- 
dimensional flows. The primary emphasis in this in- 
vestigation to date has been placed on the flow about 
bodies of revolution. A brief account of our progress is 
the subject of this paper. 

AT ZERO ANGLE OF ATTACK 


BODIES OF REVOLUTION 


The hypersonic similarity law for potential flow about 
slender nonlifting bodies of revolution was first de- 
veloped by Tsien.' He found that, through the intro- 

Presented at the Hypersonic Aerodynamics Session, Nine 
teenth Annual Meeting, I.A.S., New York, January 29—February 
1, 1951. 

* Acronautical Research Scientist. 

t Aeronautical Research Scientist. 


duction of nondimensiona!l coordinates and a nondi 
mensional perturbation potential into the equation of 
motion and the boundary conditions, the flow around 
bodies of a given thickness distribution was a function of 
a single parameter. This parameter, known as the 
hypersonic similarity parameter, is the ratio of the free- 
stream Mach 
The similarity law for nonlifting bodies of revolution is 


It should be noted 


Number to the body fineness ratio. 


summarized at the top of Fig. 1. 
that since the pressure ratio is a function of the similar- 
ity parameter, A,, the drag parameter, MJ/,? C),* is 
also a function of A; only. An example of conditions 
where the similarity law would apply is shown in the 
lower half of Fig. 1. 

In view of the simplifying assumptions involved in 
deriving the similarity law —namely, irrotational flow, 
slender bodies, and Mach Numbers much greater than 
|. the applicability of this law seemed uncertain and 
therefore was investigated.” This investigation was 
carried out by comparing pressure distributions along 
several bodies having the same values of the similarity 
parameter and the same thickness distribution. 

The most extensive comparisons were made for 
ogive cylinders. Aithough ogives do not have exactly 
the same thickness distribution, the variation was found 
to be negligible for the cases considered. Pressure 
distributions were obtained for these bodies using the 
method of characteristics for Mach Numbers and fine- 


* Force and moment coefficients in this paper are referred to 
the maximum cross-sectional area and overall length of the body 
These coefficients do not consider the effects of pressure forces 


on the base of the body 





Similarity Law 


P/Do = f(K)) 
MZCp=A(K}) where — Kk, = i 
Example si 
ti i——- 
Mog NN —_ =< 
(Vag = “ol Yadp 
Mog r Mp Mop 
Kig . Kip 
Fic. 1 Hypersonic similarity law for slender nonlifting bodies 


of revolution, 
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Fic. 2. Pressure distributions along two ogive cylinders at zero 
angle of attack for A = 0.5. 
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Fic. 38. Pressure distributions along two ogive cylinders at zero 


angle of attack for K = 2.0. 


ness ratios varying from 1.5 to 12 in combinations to 
give values of the similarity parameter from 0.5 to 2.0. 
Comparisons were made for cases in which the pressure 
determined both and with- 
that is, with and without 


distributions were with 
out consideration of rotation 
entropy gradients in the flow. 

A sample of the results obtained is shown in Fig. 2, 
which gives the pressure distributions obtained for two 
ogive cylinders having the same value of A, for the 
The pressure coefficient (p/po) — 1 is 
The value of K, 


nose section. 
plotted against per cent nose length. 
here is 0.5 with Mach Numbers, J/o, of 3 and 6, and 
nose fineness ratios, 1/d, of 6 and 12, respectively. 
The curve has been faired through the points for both 
(More points were used in fairing the curve 
It is apparent that within 


cases. 
than could be shown here.) 
the accuracy of the solution the pressure distributions 
for these two cases are identical. The similarity in 
pressure distribution holds on the cylindrical afterbody 
as well as on the nose section. For the case shown here, 
the change in the pressure distribution due to rotation 
is of the same order of magnitude as the scatter of the 


points for a given solution and, consequently, could not 
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be shown. Fig. 3 gives the pressure distribution for 
ogive cylinders for A, = 2; one case is for a Mach 
Number of 6 and a nose fineness ratio of 3, while the 
other is for a Mach Number of 12 with a fineness ratio 
of 6. 
tained by neglecting the effects of rotation, and the 
dashed line is faired through the points obtained by 


The solid line is faired through the points ob. 


including the effects of rotation. Again the pressure 
for the irrotational 
As predicted by Hayes,’ the pressure dis- 
are the 
It is noteworthy that the effect of rotation 


distributions case are nearly 
identical. 
tributions also same when rotation is con- 
sidered. 
is sizable for the case of A, = 2, whereas it was negligible 
for the case of A, = 0.5. 

These and other solutions obtained in the investiga- 
tion demonstrated that the similarity law is applicable 
to nonlifting bodies within engineering accuracy over a 
wide range of Mach Numbers and fineness ratios, in 
More 


definite limits on the range of applicability were needed, 


spite of the assumptions made in the derivation. 


To obtain these limits for ogives, the conditions at the 
nose, where the body slope was greatest, were con- 
sidered to give the most severe test of similarity. By 
following this assumption, the pressure ratios for cones 
with vertex angles corresponding to the nose angles of 
the ogives were obtained from Kopal’s tables‘ and were 
plotted against the similarity parameter for the bodies 
A curve was faired through the 
values obtained The 
similarity law was considered to be applicable for all 


under consideration. 


for the most slender bodies. 
cases where the pressure ratios were within 5 per cent 
of the faired curve. Additional solutions by the method 
of characteristics were made to substantiate the range 
of applicability as determined by the foregoing pro- 
cedure. Fig. 4 shows the range of applicability of the 
hypersonic similarity law as determined for ogives by 
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SIMILARITY LAWS FOR 
Similarity Low 
P/ Po " f (K),Kp) 
M 
Mz Cp = A( Ki, Ko) K+ we 
where 
MoG_ = B( Kj. Ka) Kp *a(1/d) 
MoCny al C( Kj, Ko) 
Example 
aa Se 
i, 
\/d = 1/2 (1/d)b . 
'e Kig * Kip 
since Mog *!/2 Mob therefore 
Ko,” Kap 
aq = 2 ab 


Hypersonic similarity law for slender inclined bodies of 
revolution. 


Fic. 5. 


this method. Within the unshaded area the hypersonic 
similarity law should be applicable within 5 per cent. 


BODIES OF REVOLUTION AT ANGLE OF ATTACK 


Potential Flow Solution 


Since the hypersonic similarity law proved to be ac- 
curate for nonlifting bodies of revolution, it appeared 
desirable to determine the form of the law when the 
angle of attack was not zero. For this analysis, the 
assumptions made were consistent with those of the 
previous investigation—that is, slender bodies, small 
angles of attack, large Mach Numbers, and potential 
flow. The equation of motion and the boundary con- 
ditions were established to the order of the approxima- 
tions involved in these assumptions and were then re- 
duced to nondimensional form. In this form, flow 
similarity was indicated for related bodies if Ay, the 
ratio of the free-steam Mach Number to the fineness 
ratio, and AK», the product of the angle of attack and the 
fineness ratio, remained constant. 

This similarity law is summarized in Fig. 5. The 
pressure ratio for inclined bodies is a function of the 
parameter, A», as well as A;. It follows then that the 
drag, lift, and moment parameters are also functions of 
K, and K.. The two cases illustrated on the lower half 
of Fig. 5 satisfy the requirements for similarity; hence, 
the pressure distributions, in terms of p/p, and the 
and moment for the two bodies 


force parameters 


should be identical. 


Viscous Considerations 


Following the suggestion of Allen,® it is possible to 
include certain viscous corrections in the aerodynamic 
coefficients of slender inclined bodies of revolution. 
For example, consider a body of revolution at angle of 
attack. The axial component of the velocity gives 
rise to surface shear, while the transverse component 
results not only in surface shear but usually in the 
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separation of the flow. This latter phenomenon has a 
pronounced effect on the lift, drag, and moment char- 
acteristics of an inclined body. It is this phenomenon 
that was considered.* 

If a body of revolution is sufficiently slender, the 
viscous flow across it may be approximately represented 
as that over a circular cylinder of infinite length. 
With this approximation, expressions were obtained for 
the aerodynamic ccefficients resulting from the viscous 
transverse force. Upon transformation into the non- 
dimensional coordinates previously employed, these 
expressions become functions of the same parameters 
obtained for potential flow. As a result, the similarity 
parameters A, and AK» apply to the viscous cross flow 
as well as the potential flow. The additional restric- 
tion, dictated by the viscous considerations, is that the 
Reynolds Number of the transverse flow must remain 


constant. 


Experimental Results 


As yet we have not been able to make a compre- 
hensive evaluation of the similarity law for inclined 
bodies of revolution. However, a limited experimental 
check of the law has been conducted in the Ames 10- 
by 14-in. supersonic wind tunnel. It will be recalled 
that the hypersonic similarity law predicts that the 
ratio of the local pressure to the free-stream static pres- 
sure at corresponding points on related bodies is a 
function of the similarity parameters only. A simple 
experimental check would then involve measuring the 
pressures at corresponding points on two cones of 
different fineness ratios and comparing the results on the 
This 


type of check has been made and is illustrated in Fig. 6. 


basis of the similarity parameters A, and Ko». 


* R. T. Jones, in N.A.C.A. Report No. 884, 1947, demonstrated 
that the axial and transverse components of the laminar bound- 
ary-layer flow could be considered separately for the case of 
infinitely long cylinders at angles of yaw. 
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3.0 and 4.9, 
were tested at Mach 4.46, re- 
spectively, and at angles of attack up to 5°. This 
combination gave a value of A, of 0.91 for both bodies. 


Two cones, having fineness ratios of 
Numbers of 2.75 and 


Pressures were measured at five points around the sur- 
face in a transverse plane at one-half the body length. 
In the lower half of Fig. 6, the pressure coefficients for 
corresponding points are plotted as a function of Ko. 
\greement with the hypersonic similarity law can be 
observed here in that the values of pressure coefficient 
for corresponding points on both bodies lie essentially 
along the same curve. A discrepancy at the higher 
values of K, can be noted on the low-pressure side of 
the cones. This is believed to be a result of dissimilar 
separation characteristics caused by the fact that we 
were unable to maintain the same cross-flow Reynolds 
Number at corresponding values of A». Separation 
phenomenon should be almost identical at the same 
cross-flow Reynolds Number. 


THREE-DIMENSIONAL FLOW 


As a result of the success in applying the similarity 
law to the specific cases previously discussed, it ap- 
peared desirable to determine the form of the law for 
three-dimensional shapes of practical value. In this 
regard, it has been indicated by Hayes* that similarity 
will be observed in hypersonic flow about slender three- 
dimensional shapes. His analysis included consideration 
of the effects of shock waves and rotation; however, it 
precluded considerations of the form of the law in 
terms of customary aerodynamic parameters. It was 
therefore undertaken to determine the law as it would 
apply to familiar three-dimensional configurations." 
The assumptions made were the same as those pre- 
viously employed with one distinct difference: The 
assumption of slenderness was applied only to the 
thickness dimension. On this basis, the boundary 
conditions and equation of motion were transformed 





Similarity Law 
p/P,* f(Kj,Ko,K3) 


M 
M3 Cp * A( Kj, Ko,K3) kK? 42 
Mo CL = B( Kj, Ko, Kz) where Ko = a(l/g) 
Mo Cy = C( Ky» Ko, Ka) Kz = 97/b 


=xample 


Mog 





(17g), "2 (Yale 


Ki, * K 
Mog =l/2 Mop le 'b 
since ag = 2a therefore Ko,* Ka, 
6 * 
(Yb), = (VYbp K3q* Kay, 


lic. 7 Hypersonic similarity law for wing-body combinations 


AERONAUTICAL 





SCIENCES—AUGUST, 109051 

into the nondimensional system of coordinates, and 
flow similarity for related bodies was then observed, 
The hypersonic similarity law for three-dimensional 
flow, together with an example of how it may be ap. 
plied to a wing-body combination, is shown in Fig, 7, 
In this case, the pressure ratio is a function not only 
of A, and A, but also A;, the ratio of the body thickness 
to the span. The force and moment parameters are 
also a function of these three similarity parameters, 
At present, it appears impractical to further generalize 
the law to include viscous effects. 

These expressions for three-dimensional similarity 
apply, as would be expected, to particular cases in 
addition to the general wing-body combination already 
illustrated. For example, consider the following cases 

(1) Two-dimensional flow. the 
A3, is zero, and therefore the pressure ratio and force 


Here, parameter, 
and moment parameters are functions of only A, and 
A». This result is identical to that developed by Tsien 
for two-dimensional flow. 

(2) Slender body of revolution.—In 
parameter, A;, has the constant value of | and is thus 
The resulting 


this case, the 


eliminated as a similarity parameter. 
expressions due to this consideration are identical to 
those for the inclined body of revolution previously 
discussed in this paper. For the zero angle-of-attack 
case, the value of the similarity parameter, A», would, 
of course, be zero. The similarity law resulting from 
this consideration (noting the fact that Ay is constant 
would be the axially symmetric case also derived by 
Tsien. ! 

(3) Finite span wings..—This is a three-dimensional 
case involving all the parameters. It should be em- 
phasized, however, that the parameter, A;, appeared to 
the second power throughout the equations of motion. 
It may therefore be concluded that only for wings where 
the aspect ratio is extremely small will the characteris- 
tics differ appreciably from the two-dimensional case. 
(It has recently been brought to our attention that this 
rule was also developed by Professor Tsien and, al 
though unpublished, has been presented in the form of 


lecture notes. ) 


CONCLUDING REMARKS 


A final appraisal of the status of our investigation of 
the hypersonic similarity law may be summarized in 
three main points. First, it has been shown that in 
spite of the assumptions made in its derivation, the 
hypersonic similarity law is applicable to bodies of 
revolution at zero angle of attack for wide ranges of 
Mach Number limited 
experimental data presented here indicate that the hy- 


and fineness ratio. Second, 
personic similarity law also applies for bodies of revolu 
tion at small angles of attack; however, additional data 
must be obtained before the range of applicability can 
Third, the of the hypersonic 


be established. form 


(Continued on page 568) 
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A Superposition Method for Calculating the 
Aeroelastic Behavior of Swept Wings’ 


R. B. BROWN,? K. F. HOLTBY,? anp H. C. MARTIN** 


Boeing Airplane Company 


SUMMARY 


This paper presents the development and application of a 
method for solving aeroelastic problems on sweptback wings 
primarily suitable for use by aerodynamicists and structural de 
sign engineers. This method, called the ‘‘superposition method,”’ 
separates the aerodynamic and structural calculations, making it 
possible to use either theoretical or experimental determinations 
of spanwise loading and to account for induction effects. The 
method has been developed principally for use in design calcu 
lations where simplicity and rapidity of solution are primary con 
siderations. 

First, certain selected initial twist distributions are applied 
to the wing. For these, the spanwise loading is determined and 
the elastic deflections are computed. Any arbitrary loading 
may then be calculated by a superposition of these initial effects, 
plus the requirement of equilibrium between elastic and acro 
dynamic forces. 

The superposition method can be used to solve the usual acro 
dynamic problems arising in longitudinal and lateral stability 
and control and to provide structural design criteria. Applica 
tion of the method to these problems is discussed in the second 
part of this paper. In applying the method, once certain pa 
rameters have been calculated, the solution toa large number of 
specific problems follows at once by application of the super 


position principle 


INTRODUCTION 


| yeas LOADS ON THE WING which cause lift, 
rolling moment, etc., are dependent on the span- 
wise distribution of angle of attack. For the nonrigid 
or flexible wing, this angle-of-attack distribution will 
be a function of wing deflection. As a result, aero- 
dynamic forces and structural deflections cannot be 
considered as being independent of each other. 

Examples of how structural deflections alter typical 
aerodynamic parameters are shown in Figs. | and 2. 
These parameters come from the solution for the 
equilibrium state between aerodynamic and _ elastic 
forces. Determination of this equilibrium state con- 
stitutes the solution to the flexible wing problem. 

In determining the equilibrium state, all deflections 
that affect the wing angle of attack must be considered. 
On straight wings, only torsional deflections cause 
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vised and received February 15, 1951 
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t Aerodynamics Unit, Boeing Airplane Company, 
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** Associate Protessor, University of Washington 


and Unwersity of Washington 


changes in angle of attack along the span. On swept 
wings, both torsional and bending deflections cause 
changes in the stream angle of attack. 

Aerodynamic loads are a function of dynamic pres 
sure, g, and section lift coefficient, C,. The lift coeffi 
cient is dependent on the wing plan form, angle-of 
attack distribution, and the two-dimensional slope 
of the lift curve, m. For a given plan form and angle 
of-attack distribution, structural deflections are directly 
proportional to gm. Accordingly, gm is a basic param 
eter that appears throughout this analysis. The vari 
ation of m with Mach Number can be taken into ac 
count either by use of wind-tunnel data or the usual 


Glauert relationship 
m = (m)ya=o (1/V 1 — M?* cos? A) 


The three-dimensional slope of the lift curve for the rigid 
wing can be used, since two- and three-dimensional 
values vary in the same manner with Mach Number. 
In this case span-load distributions must be based on a 
unit value of the three-dimensional slope for the rigid 
wing. 

Methods of solving the component aerodynamic and 
structural problems are prerequisite to establishing a 
solution for equilibrium conditions. For a given angle 
of-attack distribution, it must be possible to determine 
the rigid-wing span-load distribution. Also, for any 
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load distribution, it must be possible to calculate the 
elastic deflection of the wing. Various methods of 
handling these problems are discussed in succeeding 
pages of this paper. However, it should be noted that 
use of the superposition method is independent of the 
manner in which the individual aerodynamic and struc- 
tural problems are solved. 

Two of the earlier treatments of this problem’: ? pre- 
sent solutions based on methods quite different from 
each other and from that proposed in this paper. In 
these previous approaches to the problem, the spanwise 
load calculation is carried out for the elastic wing 
structure. Thus the span-load problem and the wing- 
deflection problem both enter more or less simultane- 
ously into the formulation and development of these 
methods. 


tion method is that the span-load problem is carried out 


The significant difference in the superposi- 


independent of the elastic properties of the wing. 
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Fic. 3. Deflection of typical sweptback wing due to lifting 
loading. 
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PART (I) DEVELOPMENT OF THE SUPERPOSITION 
METHOD 


Equilibrium Condition 


As already pointed out, the elastic wing subjected to 
aerodynamic loading will deflect until equilibrium js 
established between external aerodynamic and internal 
elastic forces. Aerodynamic loadings are caused by 
angle-of-attack distributions; also, angle-of-attack dis- 
tributions result from structural deflections. There- 
fore, the relationships upon which the superposition 
method is based are expressed in terms of angle of at- 
tack. 

The rigid wing is considered as being subjected to 
given initial conditions of angle-of-attack distribution 
a,;and qm. The aerodynamic forces thus initiated will 
deflect an elastic wing by a spanwise distribution a, 
until equilibrium is attained at a final angle-of-attack 
distribution ay. These angles are shown in Fig. 3 for 
the wing with sweepback. 

The angles shown in Fig. 3 are defined as follows: 


a, = angle-of-attack distribution initially applied 
to the wing; hence, a; refers to the rigid or 
undeflected wing 

a, = elastic angle-of-attack distribution due _ to 


wing deflection (on a swept wing both bend- 
ing and torsion contribute to a,) 
apy = angle-of-attack distribution at which equilib- 
rium is established—referred to as the final 
angle-of-attack distribution 


These angles are functions of the spanwise coordinate 
n and are measured in the direction parallel to the plane 
of symmetry of the airplane. 

The aerodynamic forces due to a, will be termed the 
rigid-wing loading. The wing will deflect causing an 
angle-of-attack distribution, a, and the loading cor- 
responding to a, will be called the “elastic loading.” 
At equilibrium, ay, the aerodynamic forces will be re- 
ferred to as either the equilibrium or the final loading. 
At every station along the span, the three angles are 
related by 


Qe = ar + az (1) 


This relationship is basic to any consideration of the 
effect of structural deflections on aerodynamic loading 
and, hence, is basic in the development of the super- 
position method. 

Another important fact is implicit in Eq. (1). Sup- 
pose the equilibrium angle of attack ais known. Then 
the equilibrium loading may be found, and, since this 
is the loading that actually holds the wing in its de- 
flected position, a, may then be calculated. Knowing 
ay and az, Eq. (1) gives a;. Hence, starting with a,, 
a complete solution is directly obtainable. 

On the other hand, starting with a,;, no such direct 
determination of the quantities in Eq. (1) is possible. 
Given a;, neither ag, nor arg can be found by direct 
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means, since they are mutually dependent on each 


other. 


Superposition of Preliminary Cases 


Superposition will be based on several solutions de- 
termined beforehand for a few selected preliminary 
cases. For these cases, ay distributions will be as- 
sumed and az and a; calculated by the procedure pre- 
viously indicated. Since calculation of a, will de- 
pend on wing stiffness and geometry, these initial 
solutions must be carried out for a specific wing.  In- 
formation regarding the wing structure must be avail- 
able if the influence of structural deflection on wing 
aerodynamic forces is to be predicted. 

Let three preliminary solutions be carried out for 
certain selected Cases 1, 2, and 3, for which a, is known. 
Then, all terms in Eq. (1) are known for each of these 
cases (i.€., @7,, Qz,, @r, for Case 1, etc.). The exact 
form for the selected values a;z,, ay,, and ay, are not 
important at this point. 

Now assume the three known initial values a;,, @,,, 
and a,;, to be multiplied by constants A, B, and C, re- 
spectively, and superimposed as shown in Fig. 4. 

The superposition of Aa;,, Ba;,, and Ca,, will yield 
which will be 


another angle-of-attack distribution, 


referred toasa,;. Or 
Aa,, 4. Ba, + Ca, = dy 


Referring to Eq. (1) and 


Aar,, = Aa, + Aas, (2a) 
Bar, = Ba,, + Bap, (2b) 
Car, = Ca;, + Caz, (2c) 


accordingly, elastic a, and final ay deflections can be 


superimposed as follows: 


A az, + Bar, + Cap, = ax 
A ar, + Bay, + Car, = ay 


Hence, by superposition of the three basic cases, a solu- 
tion a, to an arbitrary case a, has been obtained. 

The constants A, B, and C must be determined if a 
specific a; is to be duplicated. It will be the function 
of the analytical development of the method to show 
how these may be calculated. At present, it is of in- 
terest to observe that in Eqs. (2a), (2b), and (2c) all 
terms with the same numerical subscript contain the 
same constant multiplier. This is necessary if Hooke’s 
Law is to hold and if aerodynamic loads are to be linearly 
related to the angle of attack. Thus consider a,, 
a, and a,, all of which are known. 
doubled. The loading will then also be doubled. By 
Hooke’'s Law, the elastic deflection a,, will likewise be 
Finally, by Eq. (1) 


Suppose ap, is 


doubled under these conditions. 


it is found that a;, also doubles. Hence, under these 


conditions of linearity, any factor such as A, B, and C 


above must appear on each of the three appropriate 


angles. 
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Selected Preliminary Cases 


The three preliminary cases that are basic to the 
superposition method as it will be developed in this 
paper are shown in Fig. 5. These now become the pre- 
liminary Cases 1, 2, and 3, respectively, of the previous 
section. Subscript | will then always refer to the linear 
final angle-of-attack distribution, etc. 

It will be noted from Fig. 5 that a, at the tip has been 
assumed as | rad. for all three cases. This is an ar- 
bitrary but convenient assumption. 

Having chosen values of ay for the three preliminary 


cases (ay = 7, n*, n*), span-load distributions are es 
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tablished, and ay and a, corresponding to each case are 
calculated and plotted against 7 as shown in Fig. 6. 
These plots can be superimposed as indicated by Fig. 4. 


Further Examples of Superposition Procedure 


The superposition carried out in Fig. 4 was for a prob- 
lem in which the initial angle-of-attack distribution a, 
could be matched by a superposition of the three pre- 
liminary cases (a@,,, @;,, and @;,). For problems in which 
this is not possible, the use of an “auxiliary diagram” 
is introduced. The auxiliary diagram makes the super- 
position method capable of handling problems in which 
a, is completely arbitrary. The following examples 
cover cases that illustrate the use of the auxiliary dia 
gram: 

Example (1) Constant Angle of Attack. 
to be solved is given by a, on the right-hand side of Fig. 


The problem 


7. Since a, + 0 at the root, it cannot be matched by a 
simple superposition of a,'s corresponding to linear, 


quadratic, and cubic a's as was possible in Fig. 4. 
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To overcome this difficulty, another case, repre- 
sented by the auxiliary diagram, is superimposed, to- 
gether with the basic linear, quadratic, and cubic cases 
as illustrated in Fig. 7. 

The essential feature of the auxiliary diagram is that 
it is constructed from the given problem on the basis of 
putting 

iy = oly (4 
Since ay’ is known, ay’ and a,’ may be calculated pre- 
cisely in the same manner as already carried out for the 
linear, quadratic, and cubic cases. It should be noted 
that this permits the calculation of the spanwise load 
ing as though the wing were rigid. Also, for all arbi 
trary problems, the four cases shown in Fig. 7 are suf 
ficient to yield an adequate solution. However, the 
auxiliary diagram must be constructed on the basis o/ 
Ey. (3) and, hence, will be different according to the 
specific problem. Thus, a span-load problem and a 
solution for the case represented by the auxiliary dia 
gram must be carried out before the superposition 
method can be applied. This is comparatively simple, 
since the equilibrium condition is known for the auxil 
iary diagram because of its construction on the basis oi 
Eq. (3). 

Superposition now gives 


Aa,;, + Ba,, + Ca, +a; =a, (ta 
Aap, + Bap, + Car, + ar’ = ar (4b 
Aapz, + Bay, + Cap, + ap’ = ay (4e 
where ay, = 7, ar, = 7°, and ay, = 7° (see Fig. 3). 


It should be noted that, for this type of problem, An + 
Bn? + Cn* also equals ay. 

The difficulty with the nonzero value of a; at the root 
has now been overcome by introducing the auxiliary 
diagram. This can be seen by examining Fig. 7. 

Example (2) Step Angle of Attack (Atleron, Spoilers, 
Flaps). The solution of the step angle-of-attack dis 
tribution is the same as that for the constant angle-of 
attack distribution of Example (1). Fig. 8 illustrates 
the combination of a,’s which can duplicate the step 
a, The solution ay for any type of a; may be made 
if it is possible to establish the rigid-wing-load distribu 
tions for a;. 

It is to be emphasized that any other arbitrary prob 
lem may be similarly solved with the aid of the three 
basic cases (linear, quadratic, and cubic) plus_ the 


auxiliary diagram. 


Function of Auxiliary Diagram 
The true function of the auxiliary diagram can best 

be seen by writing Eq. (4a) as 
Aa,, + Ba,, + Ca;, = a; — a, (0 


Reference to the auxiliary diagram of either Fig. 7 or 
S [and keeping basic Eq. (1) in mind | indicates that 
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Therefore, the summation of initial values as given 
by Eq. (4a) is actually made to approximate az’, 
rather than a, as a result of introducing the auxiliary 
diagram. 

This is significant and 
Physically, the wing structure is 


makes the superposition 
method possible. 
such that under any applied loading it will deflect in a 
Thus, a,” will, in every case, be smooth 
Consequently, a summation as 


smooth curve. 
uid zero at the origin. 
implied by Eq. (4a) becomes entirely possible and rea- 
sonable 

Another point is worthy of mention. Other diagrams 
than those based explicitly on Eq. (3) could serve as the 
auxiliary diagram. However, there is normally no 
advantage in deviating from the procedure given above ; 
in fact, Eq. (3) gives an exceptionally convenient 
means for determining a satisfactory auxiliary diagram. 


Analytical Expressions for «1, a1,, «1 


It has already been stated that a,,, , ete., are 
known as curves plotted against spanwise coordinate 
n. These elastic angle-of-attack distributions will now 


be approximated by cubic equations in 7 as follows: 


Aap, = Gm\(ayn + hyn? + cyn*) (7a) 
ap, = gm(am + ben® + Cm") (7b) 
Ap, = gm(asn + bsn? + ¢3n°) (7c) 


where ay, ..., ete., increase with dynamic pressure q 
and slope of lift curve m. 
select three reasonable spanwise values of 7 and read 
from the plot of 


To determine a, 6, and ¢, 


off the corresponding values of a, 
ap, VS. n. Substituting these into Eq. (7a) will give 
three equations from which to determine a), 4, and 
¢. The other constants can also be found by similar 
procedure. 

Since ap, @p, and ay, are linear, quadratic, and 
cubic, respectively, in » (see Fig. 3), substitution in 


Ky. (1) gives 


a) n ar, = (1 ayqm)n biqmn 
cigmn® (Sa) 
a n” — ap, = —asqmn + (1 bogm )n? 
cogmn® (Sb) 
a n Ap, = —a3;qgmn bxgmn? + 
(1 — ¢3gm)n" (Se) 


All coefficients of » in the previous equations are 
known. Also, gm is now seen to appear as a funda- 


mental parameter in the problem. 


Solution for Superposition Coefficients 1, 8b, and C 


Now consider again Eqs. (5) and (6). Actually, 
a, 1s known as a plot vs. 7. Hence, it can be repre- 
sented also by a cubic in 9 or 

ap’ = gm(an + by? + cn? (9) 
where a, 6, and ¢ are known by an analysis similar to 
that following Eqs. (7 But by Eqs. (5) and (6) it 1s 


also seen that 


LOR OF 


SWEPT WINGS 


ap’ - as=— a," = Aa, + Ba, + Ca, 
By Eggs. (8), this last result is also seen to be a cubic in 
n. Hence, 


qm(an + bn? + cn*) — (Aa, + Ba;, + Ca;,) = 0 
(10) 


since the two cubics in 7 can be made to coincide iden- 


tically. Substituting from Eqs. (8) into Eq. (10), 


lagm — A(1 — aiqgm) + Basgm + Casgm|n + 
lbgm + Abgm — BU — bogm) + Chsqm|n? + 
|cgm + Acqm + Boogm — CU — exqm)|n§ =O (11) 


This last equation can be assured for all values of 7 


only if the coefficients of y are made to vanish. This re- 
quires 
A(l — aigm) — Baszgm — Ca;sgm = aqm_ (12a) 
—Abqm + BU — bam) — Chsqm = bqm (12b) 


—Acgqm — Bogm + C(1 — exqm) = cqm_ (12¢) 
Out of Eqs. (12) a solution can be made for A, B, and 
C. Then, from Eq. (4b), it is seen that the desired 


solution for the equilibrium a, becomes 


apy = An + By? + Cr* + (13) 


in which Eq. (3) has been used to put ap’ = a;.  Like- 
wise, ay, is known from Eq. (4c), since ay’ has been pre- 
viously determined. 

It should be noted that A, B, and C are the only un- 
knowns in Eq. (12). A specific value for the parameter 
qm must be assigned. By assuming a series of values 
for gm, the superposition coefficients A, B, and C can 
be established for the desired range of gm. 

Application of the method to a wide variety of prob 
lems has justified the use of a cubic in 9 tor approxi 
mating a,’. A higher degree curve might be used, but 
this increases the numerical work considerably. In 
certain cases trial and error may be necessary to find 
the best form of the cubic used to represent ay’. How- 
ever, a solution once found can always be checked. 
Thus, having determined a,» as from Eq. (13), it is pos 
sible to calculate the corresponding loading. Having 
found the loading, the deflections may be calculated. 
These must check the value given by Eq. (4c) using 
A, B,and C. 

Part (II) will indicate how the solution for a number 
of particular quantities follows at once when A, B, and 
C have been determined for a given problem. 
SUPERPOSITION METHOD 


ParT (II) APPLICATION OF 


As mentioned in the preceding discussion, the super- 
position method provides the designer with a rapid 
means of evaluating the effect of structural flexibility 
on aerodynamic characteristics and structural design 
criteria. It has been shown that the superposition 
method provides a means of solving for equilibrium 
conditions involving aerodynamic loads and structural 


deflections. A variety of problems requires solutions 
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Fic. 9. Lift coefficient variation 


for such conditions. Typical problems for which 
answers are frequently desired include longitudinal and 
lateral stability and control characteristics and struc- 
tural design information on wing loadings, balancing 


tail loads, structural deflections, gust alleviation, etc. 


Preliminary Information Required for Superposition 


Since all solutions are obtained by a superposition of 
known effects, the bulk of the work in using the method 
is involved in establishing certain basic information. 
This basic information is used in solving for the super- 
position coefficients (A, B, and C). Once these con- 
stants have been obtained for any given type of loading, 
a large variety of answers related to that loading be 
comes immediately available. 

As outlined in Part (I), solution for the superposition 
coefficients A, B, and C is based on combining certain 
cubic equations that represent various angle-of-attack 
distributions. The steps involved in obtaining these 
equations are as follows: 

(1) Determine the rigid-wing span-load distributions 
corresponding to linear, quadratic, and cubic angle-of- 
attack distributions. 

(2) Determine the rigid-wing span-load distribu- 
tions corresponding to any angle-of-attack distribu- 
tions to be considered which cannot be represented 
by a cubic equation (such as due to level flight, control 
surface deflection, etc.). 

(3) Apply these span-load distributions to the wing 
structure for a unit value of gm and determine the re 
sulting angular deflections along the span parallel to 
the plane of symmetry (without regard to equilibrium 
conditions). 

(4) Match the angle-of-attack distributions foun< 
in step (3) above with cubic equations of the form 
an + bn? + cn’. 

The rigid-wing span-load distributions can be deter 
mined by any one of several available methods (e.g., by 
Weissinger’s Method for sweptback wings) or from 
pressure data obtained on rigid wind-tunnel models for 
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such cases as level flight, yawed flight, control surface 
deflection, etc. 

Calculation of the angular deflections resulting from 
the application of these load distributions is relatively 
simple for wings of high aspect ratio. For these wings 
the selection of a straight elastic axis extending from the 
wing tip to the body intersection has been found to be 
reasonable. The structural problem is then reduced 
to that of a tapered beam subjected to a distributed 
load. Calculation of deflections for this 


case involves the conventional integration of bending 


structural 


and torsional moment distributions used in conjunc 
tion with stiffness values which are known as a function 
of span. The assumption of a straight elastic axis be 
comes less valid as aspect ratio is decreased. For such 
wings, a more complicated means of determining the 
structural deflections must be used. 

Since the fundamental mechanics of the superposi 
tion method are independent of the manner in which 
span-load distributions are obtained and applied to the 
wing structure, these calculations will not be discussed 
in further detail in this paper. 

Determination of the cubic coefficients a, 6, and c is 
accomplished by solving three simultaneous equations 
using three points chosen along the semispan; in this 
way it is possible to represent the angle-of-attack dis- 
tribution, computed in step (3) above, by a cubic in 
n. Choice of points to match is dependent on the shape 
of the elastic deflection curve for the particular prob 
lem being considered. 

The information listed in steps (1) to (4) is prerequi- 
site to application of the superposition method. Cer 
tain additional preliminary data is usually required, de 
pending on the particular problem for which a solution 
is being sought. For example, if one of the answers 
desired is lift coefficient C,, then C, should be deter 
mined for each of the loadings listed in steps (1) and 
(2) above. Similarly, if bending or torsional moments 
are desired as a function of span, they should also be 
computed for each of the various types of loading 
Many of these quantities are automatically obtained 
as by-products in solving for the basic information re- 
quired in steps (1) to (3). 

Table I shows typical preliminary information for 
linear, quadratic, and cubic angle-of-attack distribu- 
tions and various other types of loading which are fre 
quently analyzed when carrying out an actual design 


Solution for Level Flight Conditions 


Having determined the preliminary data discussed 
above, solution for equilibrium conditions by super- 
position is now possible for a variety of specific prob 


lems. In order to illustrate use of the method, the solu- 


— 


tion for level flight conditions will be indicated as an | 


example. The following factors cause structural de 
flections that affect the level flight lift coefficient: (1 


wing angle of attack; 


(2) built-in twist or camber dis- | 
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AEROELASTIC 


tribution; (3) basic wing pitching moment; and (4) 
wing dead weight. 
The lift coefficient for an elastic airplane can be ex- 


pressed as follows: 


BW = Cy = (Cradnce + (Cirdert Con + Cig (14) 
gS 
where 
ee lift coefficient due to wing angle of attack 
ht increment in C, due to built-in twist or 
camber distribution 
Cow increment in C, resulting from deflections 
caused by the basic wing pitching 
moment 
C.. increment in C, resulting from deflections 


caused by wing dead weight 


Each of the terms shown above can be solved for as a 
separate equilibrium condition, and the results can be 
added to give the net C,. Basie wing pitching mo- 
ment, built-in twist, and wing dead weight are usually 
Solution for 


secondary effects in the above equation. 


the lift coefficient due to wing angle of attack generally 
gives a satisfactory first approximation to the complete 
result. Fig. 9 shows a typical variation of lift coeffi 
cient due to wing angle of attack [i.e., (C,.),.,] in terms 
of the value of lift coefficient, C,,, for a rigid sweptback 
wing. 

Table II indicates the steps involved in obtaining 
solutions for each of the terms. Also shown are quan 
tities to plot which will result in single curves to be used 
in solving for any given flight condition. As indicated, 
all quantities are a function of gm. 

Using values from these curves, Eq. (14) becomes 
nil eat ty (¢ Lave m { see EN acd ol 


Oe m nm 


nm 
qs 


For any given flight condition, it is assumed that the 
following information is known: 


n = load factor = L/W 
W = airplane gross weight 
q = dynamic pressure 

S = wing area 

M = Mach Number 

h = altitude 


It is also assumed that the variaton of m with Mach 
Number is known. The only unknown in the above 
expression is, then, C,,. By solving for C,, and using 


the relation 


Cue (dC, de) rigia 


root 
the value of wing root angle of attack is found which 
represents the equilibrium attitude of the elastic air- 
plane in level flight. Cases in which the zero lift angle 
is not zero are automatically taken care of by including 
the increment in C, due to built-in twist or effective 
twist due to camber. 
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One of the most important elastic effects on swept 
wings is the shift in spanwise center of pressure caused 
by deflections. With a sweptback wing in level flight, 
as speed increases, the center of pressure moves inboard 
along the quarter chord line and, consequently, for 
ward. This shift is of primary importance in longi 
tudinal stability calculations and structural design for 
balancing tail loads. The expression for net pitching 
moment coefficient is analogous to Eq. (14) given pre 


viously. It can be written as follows: 
C, = a. + Cae + Cas + Ce. (19) 
where 
(Cyadun = the pitching moment coefficient about 


caused by 


the M.A.C. 
placing the wing at an angle of attack 


one-quarter 


(Cyr), = the increment in Cy due to built-in twist 


(Cy)aa = the increment in C,, due to the basic 
wing section pitching moment 
Cu. = the increment in C,, resulting from the 


deflection of the wing under the dead 
weight load 


As before, each of the above terms is solved for as a 
separate equilibrium condition, and the results are 
added to give the net Cy. 

Using values from the plotted curves of Table II, 


Eq. (15) becomes 





In 
fli 





1951 


AUGUST, 


SCIENCES 


AERONAUTICAL 


THE 


OF 


JOURNAL 


° 
NS 


Y~ee YF SS 


Y =m 


=. Tw ~~ FO 


mk 
the 


Ob 


Co 

in 
lat 
sol 





by 
tor 
are 
rat 





















































“D2 “Q “vo 
— 
—— | ee5 = o(f-B) “> tto“g-ow-|5 % LHD1/3M 
ok yO + 'Y4y =“%Y wd Sw Wa My Hy ‘Wainy = “Ms ao 2894) -(29- 4b) “9+ 'q“7- % 9 *% avzsa 
“D+ *OO+ “Pry="B = S7545, 2 “Ia try = “7 iy w oe > ' . aaa 
£6'D+'26%9+ O'y=‘6 sd ° vere wee xe ” 04)-* 9-(v - 4B)“ 7 9 
See aes ww Ww 
Sys yy “Oe YW “py = “Ww ~ cA a” ain 
rot "9+ Wa y= "ey “ih liattiiitided 3 wpb)" + 2D ->'y-| 9% — | IWawow 
Sy yy ty Oy yry =y — "04+ WD = “> atm = §9%5-(89- 48) 94 9 17- %> % *%b nie ae 
5p “> e~-O+ O-y="6 “D0 fraw z "Pm Wy nd Sp w w, f ’ 2 ’ ” 
fa"D+*#0"O + 9 %y="6 Doe +" "y=" |asme= FON D-wg-(v-4b)"Y| 9 9 'o 
ad +39 +q +y 
up 
a &9| 49 «(to -41B) 49 + Foi - 'o4+y-|f> | fp 
wb oe , ates LSIML 
yey sae ytin ty = my)! abate noaaratn am h9)| 4a 08 An - (89-1 '94y-| > 
90+ "09+ b+ + Da (-Y) ae © ideas 
29 Oe 4949474479 =D) | +d = SD4D - AG ('v- wh 2 '9'0 
°9 4027949 + 6417+ 46 =") ~ ial y za ) e ( 4h) ty 9g 
ey Per 7 
ae oy *® 
2 t ds, 38, 
af sind-tlafad y= 4) il . 5 €2 wb "5442p - 'o"y-| > "9 Sp v 
PW OH Wr r=) = (o- 2) 9 - WILLY 
w 
yD YD + Yi "4 = “PT"Y) is % Silla ad al *q =8q%)-(0q-4B)*94 '9 *y-| 92> "q 8B a 
97+ P+ O74 B= MEY) | ee v7 297] 2 € »  éeoear”l 
HD Be D479 = 7)" d = fb) - "B-('v- HB 29 '0 
£0") + 99 + 0"7+"B= ""9) htidnibututdiciilie i. ‘oO: ar 9 
Tw7awL WworF 
(7 a7evL Wosy euva ON/S?) LO7ed (7 a7e@vL WORX viva OWsSn) 98 @‘vy SLNWD7SITOD SINI/D/IIFO? 9N/AVO7 
SLINSFA A2v/7/KNV “22 7D 2Oz NOULNTOS POF JAWS QL SNOLLVODI |o—gno wMowy| FO FSAL 





SNOLLIGNOD LHOITS TINIT HOA NOUNTIOS 


IT F79VL 
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i. (Ca ot m+ (Cur) net m + (Cap) nce + Cum nm 
m m nm 
In this case, Cy is the unknown quantity for any given 
flight condition. 

Many other quantities associated with level flight 
can be determined by combining the proper values 
with the superposition coefficients A, B, and C. Typi- 
cal quantities for which solutions are frequently de- 
sired are shown on Table II under the heading of 
“Auxiliary Results.” 

While this paper is concerned only with wing effects, 
complete airplane analysis must include the effects of 
body and tail flexibility. Tail flexibility is analogous 
to wing flexibility. Wing, body, and tail equilibrium 
conditions are coupled and must obviously be treated 
as such for a complete solution. In most cases, how- 
ever, wing flexibility is the largest contributing factor 
in any consideration of elastic effects. 


Solution of Lateral Control Problem 


Determination of steady-state roll conditions will be 
discussed as an additional example of application of the 


superposition method. 


Steady-state roll occurs when the roll damping 
moment equals the moment produced by the control 
device. For rigid wings, the damping load is a func- 
tion of the roll angle of attack along the semispan 
(pb/2v)n. This corresponds to a linear angle-of-attack 
distribution along the semispan having a value of pb + 
2v at the tip. The helix angle of rolling flight is also 
equal to pb/2v. Consequently, pb/2v becomes the com- 
mon parameter for comparing roll characteristics. 
Equations relating damping and control moments are 


as follows: 


| a = (dC, dé) dqsb 


dC, pb 
M temsies = qsb 
d(pb/2v) 2v 


and defining 


P dC, . dC, 
L3 = -— i = 
7 d6 . d(pb 2v) 
then, since DF sacar = M samping (for steady state), 


pb/2v = (CL; C1.,)6 


Obviously, only unsymmetric loads can cause roll. 
Consequently, symmetric loadings due to basic wing 
angle of attack, dead weight, etc., do not enter into the 
lateral problem. The problem therefore reduces to 
solving for CL, and Ci, as quantities that are influenced 
by structural deflections. Typical plots of these terms 
for a swept wing with an aileron-type control system 
are shown in Fig. 10. 

Table III gives the solution for these terms as sepa- 
rate equilibrium conditions. The rigid value of as = 
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da/dé6 must be known or estimated for the particular 
control configuration under consideration. 


CONCLUDING REMARKS 


Having obtained the superposition coefficients A, B, 
and C for any given type of loading, it is always pos- 
sible to perform a check to determine whether or not 
these coefficients actually correspond to equilibrium. 
When the wing is subjected to the final equilibrium 
loading corresponding to ay, it must cause the struc- 
ture to deflect to ay. The deflection computed from 
application of this loading should then equal An + 
By? + Cr’. 


C have been computed incorrectly or the various cubic 


If a discrepancy exists, either A, B, and 


coefficients used in solving for A, B, and C do not 
match the deflection curves with sufficient accuracy. 
On the wings examined to date it has been found that 
the latter difficulty can be eliminated by careful selec 
tion of matching points. 

It should be noted that as gm increases it is necessary 
to increase the accuracy of the computations in order 
that A, B, and C can be closely checked. However, 
the closeness of the check does not indicate the accu- 
racy of A, B, and C unless the gm for which they were 
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determined is considered. This is due to the fact that, 
at high values of qm, slight deviations from the true 
equilibrium angle of attack, a», will produce large loads, 
which in turn have a big effect on elastic deflections, 
This is shown by Fig. 11. 

Application of the superposition method is not lim 
ited to solution of the problems shown as examples on 
the preceding pages; nor is the method limited to use on 
swept wings. The relative simplicity of solution makes 
this method particularly adaptable for use in airplane 


design. 
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Nonsteady Supersonic Flow About Pointed Bodies of Revolution 


Continued from page 511 


70 
I, = / S’(x — BR cosh uv) cosh u du (Al) 


cosh (x/BR 


S’(x — BReosh u) X 


“0 
IE = / 
27 cosh (x BR) 


(x — BR cosh u) cosh udu (A2) 


1) 
I; = / S(x — BR cosh u) cosh u du (AS 
J cosh “x BR 
0) 
I, = | S'(x — BR cosh u) du (A4 
cosh ‘(x /BR) 


When these integrals are differentiated with respect 
to x, it should be remembered that S(O) = S’(0) = 0 for 


pointed bodies of revolution. 
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Transient Loading of Wide Delta Airfoils at 
Supersonic Speeds 


JOHN W. MILES* 
l’nwersity of California at Los Angeles 


SUMMARY 


The indicial admittance of a pointed wing, bounded by two 
supersonic leading edges and a straight trailing edge and located 
in an otherwise uniform supersonic flow, is calculated on the 
basis of linearized theory. This result is applied to the calcula 
tion of the lift and moments on a ‘‘wide”’ delta airfoil due to a 
sudden change of angle of attack, a uniform, sharp-edged gust, 
an impulsive gust, and sharp-edged and impulsive gusts having 
utisymmetric linear velocity distributions over the span 
The response of a two-dimensional airfoil to an impulsive gust 
is also calculated. These results are given in the form of alge 


braic formulas and curves 


NOTATION 


1 = indicial admittance, Eqs. (2.1) and (2.4) 
A = two-dimensional indicial admittance, Eqs. (2.5) and 
(2.6 
Cc) = lift coefficient 
Ci, = Cx due to unit increment in angle of attack, Eqs 
(3.8) and (3.11 
CE = (C;, due to impulsive gust, Eqs. (5.4) and (5.6) 
Ci, = C, due to uniform, sharp-edged gust, Eqs. (4.9) and 
(4.11) 
Cu = pitching moment coeflicient 
Cu, = Cy due to unit increment in angle of attack, Eqs 
(3.9) and (3.12 
Cw, = Cy due toimpulsive gust, Eqs. (5.5) and (5.7 
Cy, = Cy due to uniform, sharp-edged gust, Eqs. (4.10 
and (4.12) : 
Ci = rolling moment coefficient 
Ci, = (; due to impulsive, antisymmetric gust, Eqs. (6.8 
and (6.9) 
Cy = (;)due to sharp-edged, antisymmetric gust, Eqs. (6.5 
and (6.6) 
G = Green's function for harmonic motion problem, Eq 
(1.8 
‘id = lift force 
J = pitching moment about a point ac /2 aft of the apex 
WV = rolling moment about mid-span 
S = surface area 
S = see Eq. (3.6) and Appendix 
1 = see Eq. (4.6) and Appendix 
= free-stream velocity 
= fraction of root semichord of pitching axis aft of 
apex; also sonic velocity 
= root chord of delta wing (or simply unit of length in 
more general use 
= spanwise weighting function, Eqs. (1.3) and (1.18 
R = reduced frequency, Eq. (1.9 
" = two-dimensional] lift curve slope (= 4 tan @ 
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p = gage pressure 

t = time 

ty = see Eq. (1.10) 

i = dimensionless time, Eq. (1.11 

u = see Eq. (4.3) 

My, = see Eq. (2.7) and (2.8) 

Ww = dimensionless (i.e., divided by 1’) downwash at air 
foil 

x, ¥, 2 = coordinates (Fig. 1) 

Fos = port and starboard leading edges 

a = angle of attack 

6 = Dirac delta function, Eq. (5.2) 

on = see Eqs. (3.17), (3.18), (3.19), and (4.13 

n = running y coordinate 

A = Machangle, Eq. (1.2) 

iv = semiapex angle 

é = running x coordinate 

p = air mass density 

T = (Ut/c) (1 — sin @) 

o = velocity potential 

x = see Eqs. (1.3) and (1.13) 

w = angular frequency 

1(x = Heaviside step function (= lifx > 0; =0ifx < 0 


(1) INTRODUCTION 


|» A PREVIOUS PAPER! we considered the aerodynamic 
forces on pointed wings, bounded by two supersonic 
leading edges and a straight trailing edge transverse 
to the main flow, executing harmonic motion in an other- 
wise undisturbed, supersonic flow. It was found that 
the problem could be reduced to an essentially two- 
dimensional one, which was particularly simple in the 
case where only the lift, pitching moment, and rolling 
moment were required. Explicitly, if the main flow 
(L’) is directed along the positive x axis, the nose of the 
airfoil projects on the origin, and the port and starboard 
leading edges are specified by y,(xv) and y,(x), respec 


tively, subject to the restrictions 


yp(O) = 9¥.(0) = O (1.la) 

y,’(x) > tan 6 (1.1b) 
Vp (x) < — tan 0 (1.1e) 
6 = sin“ (a/U) (1.2) 


6 is the Mach angle, and a is the sonic velocity. If 
g(x, y, 2) exp (twf) is the velocity potential, w(x, y) 
exp (twf) is the normal velocity (downwash) at the 


wing, positive down and referred to l’, and x(x) 1s a 
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weighted spanwise integral of the potential defined byt 


¥,(*) 
x(x) = UM f (x, vy, OF)f(y) dy 


(x) 
p 


(1.3) 


a solution to the (linearized) boundary-value problem 


l o 0 \? 
v% =—(U ) (1.4 
° ={ Ox ™ Ot ? ) 
o:(x, y, OF) = —Uw(x, y)1(x)1[y.(x) — y]1 ly — yp(x)] 
(1.5) 
(x, y, OF) =O, x =D (1.6) 


is given by 


x ™y, (&) 
x(x) = tan of arf dnG(x — & n)w(é ) (1.7) 
0 ¥,(&) 


? l 
G(x, y) = exp (—ikx sec? #)1(x)-— X 
T 


/ cos (kx sec 6 tan @ sin ¢)-f(x cos ¢ + y) dg (1.8) 
0 


k = (we/U) (1.9) 


where ¢ is a unit of length, all coordinates being dimen- 
sionless and referred to this length, and & is the re- 


duced frequency. 


In the present paper, we shall be interested in the 
response of this same wing for the case where the 
downwash w(x, y) at (x, y) arises abruptly at fo(x, y), 
and is thereafter constant in time, whence the boundary 
condition, Eq. (1.5), is replaced by 


— Uw(x, y)1[t — to(x, y)] X 
1(x«)1[y.(x) — v]l[y — y,(x)] 


(x, y, OF, t) = 
(1.10) 
where we now introduce the explicit dependence of 
gont. Eqs. (1.4) and (1.6) remain as stated. 

In considering this problem, we shall find it con- 
venient to introduce the dimensionless time /, defined 
by 

i = (Ut/e) (Eth) 


Moreover, if c is the wing chord, 
In the subsequent 


noting that wf = i. 
?/ measures time in chord lengths. 
sections the use of this dimensionless time will be im- 
plicit, and the bar will be dispensed with. 

Having found the velocity potential on the upper 
surface of the wing—i.e., ¢(x, vy, 0*, ?)—the gage pres- 
sure is given by 

D : 


(x, y, OT,.¢) = 
pi? ; 


’ ’ ’ ola 


defined such that 


p(x, y, 1) = —p 


t In reference 1 we introduced y(x, &), 


x 
x(x) = tan of v(x, &) dé 


AERONAUTICAL 


SCIENCES—AUGUST, 1951 
and a weighted spanwise integral of this pressure may 
be expressed in terms of x(x, ?) in accordance with 


y,(*) 
-f p(x, y, Df(y) dy = 
(242) x09 ct 
p a x(x, f) (1.13) 


(2) INDICIAL ADMITTANCE 


The response of a linear system to the Heaviside 
excitation 1(f) is designated as the indicial admittance 
of the system, whereas we may define the response of 
the same system to the harmonic excitation exp (ik?) 
as the transfer function. If these two characteristic 
functions are denoted by A(f) and G(k), respectively, 
then A is the Bromwich integral of G—viz., 


l 
Ali) = lim 


oT e—>0) 


} dk(e + tk)~' exp (tkt)G(R) 
(2.1) 


by virtue of the Fourier representation 


7 
1) = = lim / dk(e + tk)~' exp (ikt) = (2.2) 


<7 ~—>0 

These relations enable us to transform the solution 

to Eqs. (1.4), (1.5), and (1.6) given by Eqs. (1.7) and 

(1.8) to the solution to Eqs. (1.4), (1.10), and (1.6) 
Thus, if the latter solution is represented by 


x(x, t) = tan of ae f dnA X 
0 y,, (&) 


[x — & , t — to(&, n) }wlé, n) = (2.3) 


the indicial admittance is obtained by substituting 
Green's function, Eq. (1.8), in Eq. (2.1). 


trum integrations may be effected directly from Eq. 


The spec- 


(2.2), and the desired result is 


ee wz 
A(x, y, t) = 1(x) - FE | dg f(x cos ¢ + y) - 
oT ( 


}1[* — x sec? A(1 + sin 6 sin g)] + 1 X 


[{ — x sec? (1 — sin@sin ¢)]{ (24) 


We shall be interested further only in the special 


cases fo(y) = 1 and fi(y) = y, corresponding to the 
calculation of the lift and rolling moment, respec 
tively, on the spanwise strip at x. Substituting fy) = | 


in Eq. (2.4) and carrying out the ¢ integration yields 


foly) = 1: A(x, y, t) = Ao(x/t) (2.5) 


l(u)l(u, — u) + l(u — m,)1 (te — u)- 


| u — Uju2 | 
cos”! (2.6) 
T (l — m)u 


A of u ) = 


u, = | — sin é (2.7) 
uw = 1+ sind (2.8) 

The result obtained for f; = y may be written 
fly) =»: A(x, y, t) = yAo(x/t) (2.9) 





su 


W 
ch 
all 


M 


Int 
wri 


ciel 
len 


+ 


unit 


may 


side 
ince 
e of 
tkt) 
istic 
ely, 


ion 
ind 
6) 


ial 
he 


9) 














TRANSJENT 


The indicial admittance Ao may be recognized as 
that for the two-dimensional problem, which has been 
discussed adequately in the literature,’~® including, in 
particular, physical interpretations of the waves repre- 
sented by the two step functions in Eq. (2.4). More- 
over, in the important special case where fy is inde- 
pendent of y, the solution represented by Eqs. (2.3) and 
(2.4) may be cast in the form 


. 


x —— t 
x, (x, ft) = tan 6 / Ay |, e] w,(é)dé& (2.10) 


*y (x) 
w, (xX) -f wx, yy" dy, n= 
V¥,( xX) 


which reduces the problem to an equivalent two- 
dimensional problem of downwash distribution w,,(.x), 
as was to be expected from the analogous result for 


0,1 (2.11) 


harmonic time dependence. ! 
We now proceed to consider important examples 
that are covered by Eqs. (2.10) and (2.11). 


(3) SuDDEN CHANGE IN ANGLE OF ATTACK 


In the case of a sudden change in angle of attack or 
sudden control surface deflection at time zero, we have 


i(x) = 0 (3.1) 


We shall consider here only the case of a sudden (unit) 
change of incidence (without rotation) over the entire 
airfoil, whence 


wx, vy) =1 (3.2) 


Moreover, we shall assume a delta plan form, for which 


y,(x) = —y,(x) = xtanp, pw 8 (3.3) 


where y is the semiapex angle, as shown in Fig. 1. 
Because of the symmetry of the downwash distribu- 
tion |Eq. (3.2)] with respect to mid-span (y = 0), we 
shall be interested only in the lift and pitching moment 
and require only xo(x, ¢) [see Eqs. (2.10) and (2.11)]. 


Substituting w = 1 and f) = O in the latter equations 
yields 
i x—& — 
xo(x, 2) = 2 tan pw tan 6 Ay Edé (3.4) 
Introducing the change of variable (x — &) = tu, we 
write 


2 tan uw tan 6f?[(x/t)So(x/t) — Si(x/t)] (38.5) 


/ ; u"Ao(u) du 


In calculating the desired lift and moment coeffi- 


cients, we choose the root chord? (c) as our unit of 


x(x, t) = 


3.6) 


S,(z) = 


length. The surface area is given by 


S = c* tan pg (3.7) 
+ In refererce 1, the average chord (b = c/2) was chosen as the 


unit of length 
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By referring to Eq. (1.13), the lift and moment (about 
an axis adc/2 aft of the apex) coefficients can be written 


(3.Sa) 


pl? 1 "1 *y (x) 
= ( > tan x) / ax dy X 


[—2p(x, vy, O*)] (3.8b) 


7 fee 
= 4 cot yu ax + % (x,.0) (3.8¢) 
J 0 Ox 


Ot/ 
: | ks : 1 
Cur, (4) = (= .s.) M, 


U* tan ie | 
‘ (’ ‘) J dv( = «)x 
| 0 2 
/ . dy [—2p(x, y, 0*)] 


t cot ff ea 2x) (? +°) (x, t) 
= 4co : — =. x, 
eS ox or) *° 


(3.9¢) 


3.9a) 


3.9b) 


where we have used the subscript a to denote the unit 


increment in incidence. Carrying out the required 


differentiation of xo(x, ¢) in Eq. (3.4) yields 


[(0/Ox) + (0/Of) |xo(x, 4) = 2 tan pw tan 6 X 

[(x + t)So (x/t) — 2tSi(x/t)] (3.10) 
which result may be of direct interest in calculating 
chordwise loading. Substituting Eq. (3.10) in Eqs. 
(3.8c) and (3.9c) and integrating repeatedly by parts, 
we may place our results in the form 


m{(1 + 2t)So(1/t) — 
2t(2 + t)Si(1/t) + 3S2(1/8) | 


Crt) = 
(3.11) 
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@) = 6 8 
eg 
Fic. 3. ¢( Ut/c) vs. r [see Eq. (3.18) ]. 


4 . 
= —m4 (, — a} + 2(1 — ajt| X 


So(1/t) + 2tlat — 201 — a)JSi1/t) — 


Cu, (t) 


8 
t?(3a + 2t)S2(1/t) + 3 t3.$3(1 ny (3.12) 


m = 4tan 0 (3.13) 


where m is the steady-state, lift curve slope. 

The required S,(z) are evaluated in the appendix, 
and substitution of the results in Eqs. (3.11) and (3.12) 
yields 
l 
m 


; l 
Cr (t) = {(t) + cos 0 (: + 2 sin? a) x 


l ; 
M(t) + 53 — 0) f(t) (8.14) 





ae | | | 


ia 




















| 
0 2 4 
" 


) vs. r [see Eq. (3.19)]. 


ee. 2 
Cut) |a=43 = - P(1 — 2) sin® @ cos 6 H(t) + 
m 


[((1 — t)? + 2¢ sin? O]fo(t) (3.15) 


3.16) 


. A. 
= Cu lh) « ist (O-“) ao 


»? - 


3.Ld¢a) 


0O<¢t< (1 + sin @)~! 


] 1 — ¢t cos" 6 
= cos! e ) ’ 
T sin @ 


(1 + sin ¢)—' <7? < (1 — sin @) (3.17b) 

=] @> (1 — sm) (3. 17¢) 

(ft) = | (3. 1Sa) 

= (l/r) cos! [(¢ — 1) (f sin @)]} (3. 1Nb) 

= 0) (4.1Sc) 

fo(t) = O (3.19a) 
= (cos 6/r) [1 — (1 — sin @)t]’* X 

(1 + sin @)t — 1] (3.19b) 

= 0 (3.19¢) 


where the time intervals in Eqs. (3.18) and (3.19) 
correspond to those of Eq. (3.17). 

In plotting the results, it is convenient to refer them 
to the time required to reach the steady state—viz., 
(1 — sin 6)~-'. Hence, if ¢ is the real time [see Eq. 
(1.11)], we introduce the dimensionless time 
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7 = i(1 — sin 6) = (Ut/c) (1 — sin 6) (3.20) 


The functions fo, 61, ¢, (1, m)C, r and (1/m)Cy _, as 
c a 

given by Eqs. (3.14) to (3.19), are plotted as functions 
of r in Figs. 2 to 6 for various values of 1. 
(4) UnirormM Gust ENTRY 

As a second application of practical interest, we 

consider the entry of a sharp-edged uniform gust of 

unit magnitude. Such a situation is represented by 


Wx, y) = (4.1a) 


io(x, v) = x (4.1b) 


and is illustrated in Fig. 7. 

Substituting fo(&) = & and m = OQ (since we shall be 
interested only in the lift and pitching moment) in Eq. 
(2.10), we have 


‘ * a4 
2 tan uw tan 6 Ay ; 2 Edé 
“0 —_ g 


. 


xo(x, 1) = (4.2) 


To carry out the required integration, we introduce the 


change of variable 


u= (x — &)/(t — &) (4.3a) 
—& = (x — tu)/(1 — wu) (4.3b) 
t= [x — (1 — w)é]/u (4.3c) 


Because of the singularity in this transformation at 
u = | (t = x), some care must be exercised in selecting 
the limits of integration for uv. Referring to Fig. 8, 
where we have given a graphical representation of the 
transformation, we observe that the path of integra- 
tion for u, corresponding to the path 0 < & < x, de- 
pends on whether ¢ is greater or less than x. 


Ct) > 12 < 25 2 >(x/}) > u > 0 (4.4a) 
Giey <i: = F< ze > 
(x/t)<u< +a = (4,4b) 
(—eo < # < 0) 


However |see Eq. (2.6)], Au(w) vanishes for u > us and 
for negative u, and the latter path may be reduced to 
(Iv t)<u < Ue. 


Accordingly, we may rewrite Eq. (4.2) in the form 


*x/t 
2 tan pw tan O(t — x) / 
u2l(x —1t) 


(1 — u)-*Ao(u) du 


%elx, 2) = (x — tu) X 


(4.5) 


where the lower limit is 0 or #2 as the upper limit is less 


or greater than unity. Introducing the integrals 


7 ais) - | (1 — u)~"Ao(u) du (4.6) 
nol (z 1 
we write 
xo(x, f) = 2 tan pw tan O[t(t — x)72(x/t) — 
(t — x)?73(x/t)] (4.7) 


The lift and moment coefficients can now be calcu- 
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Fic. 5. Cr Ut/c), the lift coefficient due to an abrupt unit 


increment in angle of attack, vs. r 
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Fic. 6 
troidal axis, two-thirds of the root chord from the apex) coef 
ficient due to an abrupt unit increment in angle of attack, vs 
v. 


Cu (Ut c), the pitching moment (about the cen 
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Fic. 7. Wing penetrating sharp-edged gust 
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(4.3), showing (shaded area) domain of Ao(1). 
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Fic. 9. ¢3( Ut/c) vs. r [see Eq. (4.18) ]. 


lated from Eqs. (3.8) and (3.9). Carrying out the re- 


quired differentiation of xo(x, f) yields 
[(0/Ox) + (0/0f) |xo(x, t) = 

2 tan 6 tan p(t — x)T2(x/t) (4.8) 
Substituting this result in Eqs. (3.Sc) and (3.9¢) and 
integrating by parts [noting that the discontinuity in 
the derivative of 73(z) at z = 1 introduces no difficulties 


Graphical representation of the transformation, Eq. 


AUGUST, 1951 


due to the factor (¢ — x)* which arises in the partial 
integration ], we may express the results in the form 


(1/m)Cz,(t) = PSo(1/t) — Cl — 4)?T2(1/t) (4.9) 


= l 2 : 
Cu, (¢) = re Cz, (t) = 
m m >) f 


oF ue on. f* 
jes ( ) ~~ (i = "Tr ( } (4.10) 
3 t i 


where we have introduced the subscript 1 to denote the 
unit gust. 

The required integrals are evaluated in the appendix, 
and substitution of the results in Eqs. (4.9) and (4.10) 
yields 


Cr, (t) == Fo(t) oe t? cos 6¢1(t) == (1 = t)*&3(t) 


m 
(4.11) 
= 4, , 
Cu, (t) = 5 “(1 — t) cos 6¢,(t) + 
m 3 
2 2 
» Pelt) + Stl — bt) Fst) (4.12) 
» oO 
e(t) = 0 (4.13a) 
(1 — ¢t) cot 0 t sin 6 
— cosh ~! (4.13b) 
T t— 1 
= 6 (4.18¢) 


where the time intervals in Eq. (4.13) are those of 
Eq. (3.17). 
The functions ¢3, (1 m)Cz,, and (1/m Cu, are plotted 


vs. Tin Figs. 9 to 11. 


(5) IMmMpuLsiIvE Gust 


The sharp-edged uniform gust dealt with in the pre- 
ceding section may not be realized with any frequency 
under actual flying conditions (although it is, of course, 
true that a gust with any space-time history can be 
built up by a Duhamel superposition of sharp-edged 
A situation that is likely to be realized rather 
Such a 


gust is specified by the Dirac delta function—viz., 


gusts). 


more closely in practice is the impulsive gust. 


w(x, t) = 6[(Ut/c) — x] = 6(f — x) (5.1) 


/ nish as = 1, «> 6 (5.2a) 


6(z) = (d/dz)1(z) (5.2b) 


where Eqs. (5.2) define 6(z). Since w is dimensionless 
and referred to U’, it follows from Eqs. (5.1) and (5.2a) 
that the velocity time integral (at a fixed point on the 


airfoil) divided by the chord length is unity—viz., 


ae ar ; : 
/ Uw(x, t) dt = / 6(f — x) di = 1 (5.3) 
Cc . 


It follows from Eq. (5.2b) that the response of the 
airfoil to an impulsive gust is obtained from the re- 
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sponse to a unit gust after differentiation with respect 
to the dimensionless time /, which is there denoted by 
i. Carrying out the differentiation of Cy, (4) and Cu, (4), 


as given by Eqs. (4.9) and (4.10), yields 


Cr (t) = m[2tSo(1/t) + 20. — )72(1/t)] (5.4) 
Cu, (t) = m[2(a — t)tSo(1/t) — 4t2S,(1/t) + 
2(1 — t) (a — 1 — #t)T2(1/t)] (5.5) 


(The terms involving the derivatives of So, S;, and 7% 
cancel identically.) The subscript 6 denotes the im- 


pulsive gust. Substituting So, Si, and T2, we obtain 


Ci (t) = 2t cos 6 f(t) + 2¢3(t) (5.6) 
m 
l 2 
Cy, (0) = 4 (<1 ~ r) cos 66;(t) + 2tfo(t) + 
m 0 
z 
a (l — 31) &3(t) (5.7) 
o 


These results are plotted in Figs. 12 and 13, respec- 
tively. 


(6) ANTISYMMETRIC (ROLLING) GuSsTS 


The previous situations were symmetrical and, there- 
fore, did not give rise to any rolling moments. It is 
also of some interest to study the response of the wing 
to an antisymmetric excitation. As a representative 
(although highly idealized) problem, we shall consider 
the case 


w(x, y) = ycot u (6.1a) 


f(x, Vv) = x (6. 1b) 


Viewed spanwise, this distribution is shown in Fig. 7, 
while a chordwise view is shown in Fig. 14. 

Since the antisymmetric gust under consideration 
does not give rise to any lift or pitching moment, we 
shall be interested only in the rolling moment and, 
Substituting the distribution 
1 therein 


therefore, only in x,(x, /). 
(6.1) in Eqs. (2.10) and (2.11) and taking n = 


‘s “x — 
~ tan 6 tan? yw / Ay ( - 
3 J 0 t—é 


Introducing the change of variable (4.3), we 


yields 


xXi(x, /) = (6.2) 


obtain 


») 
x(x, t) = = tan’ uw tan 6(t — x) X 


7) 


» 


j (x — tu)? (1 — u)*Ao(u) du 
uel(x—t 


This result evidently can be expressed in terms of the 
T’,(x/t) [see Eq. (4.6)], as can its substantial time 


(6.3) 


derivative. 


The rolling moment coefficient is defined by 
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Ci,(Ut/c), the lift coefficient due to penetration of 
a uniform sharp-edged gust, vs. r. 


Fic. 10. 
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Fic. 11 
axis, two-thirds of the root chord from the apex) coefficient due 
to penetration of a uniform, sharp-edged gust, vs. r. 


Cy, ( Ut/c), the pitching moment (about the centroidal 


U? 1 ’ 
C4) = ¢ -S-2c tan u) M,(t) (6.4a) 
-. ig te 
= —(pl* tan p)~! / dx / y dy X 
[—2p(x, y, O*)] (6.4b) 


. 
1 


O 
= —-2 cotta | ax( > +2) xi(x, 0) (6.4¢e) 
0 g 


Substituting xo(x, /) from Eq. (6.3), carrying out the 
indicated differentiation, and integrating repeatedly 


by parts yield 
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] l ] ] l 3 l | l 
C(t) = - t4So ) (1 — t)?T» ) t(1 — t)%3 + ( ) — t)*T ( ) (6.5 
zy 3 (- +5 ) ( +; ( t)37 +, 0 t)'T, 6.5) 


Substituting the required integrals from the appendix, we have 


yn ; , 2 ‘ 
— ( ) Ci, (2) = (o(t) + ?t* cos 0¢,(t) + (1 —?t)? (1 a ‘ esc” 0¢5(t) tee 


m 
” . ” l 9 » « 
(1 — #)?] (80 + 2¢ +1) + = (1 — t)? esc? 8] &3(t) (6.6) 
This result is plotted in Fig. 15. 
It is also of interest to consider the antisymmetric impulsive gust, specified by 
w(x, y, t) = y cot pé(t — x) (6.7) 


Following the argument of Section (5), the response to this gust can be obtained by time differentiation of the 
response to the gust (6.1). Differentiating Eq. (6.5) with respect to time yields 


2 m{f.. fl ot 9 sare, ON 
C,,(t) = - 3 }#s.(‘) + 3t(1 — t) (24 — 1)7T2 (') + 2(1 — ¢t)*(4t — 1)73 (') + 3(1 — t)*74 () (6.8) 


are : ; dD l 
_ Ci,(t) = f° cos 0,(t) — 5 fC — ?t) csc” 062o(t) + (1 — 2) 3 + : (1 — t)* ese? | ¢3(t) (6.9) 
m 2 7 


This result is plotted in Fig. 16. 


(7) Two-DIMENSIONAL GusT PROBLEM ra) re) x i 

( - ) X) (x, 1) = tan 0Ay ( ) (7.3) 
For the sake of completeness and in order to find the Ox — Of t 

response to an impulsive gust, we shall consider briefly The lift and moment coefficients are given by 

the two-dimensional gust problem. In this case, we 

simply set w,(£) = 1 and fo(£) = &in Eq. (2.10), whence 


Da il il 0 Oo . 
x(x, ¢) = tan 6 f Ae (; = :) dé (7.1) = if dx (2 +2) w(x, £) (7.4b) 


= mtS)(1/t) (7.4c) 


Ci, (t) = [(eU?/2)-c]“'L (7.4a) 


Introducing the change of variable (4.3), we have 


“x/t Cy, (t) = [(pU?/2)-c-c]-'M (7.5a) 
x(x, t) = tan A(t — x) (1 — u)—#Ao(u) du (7.2) 7" | , 


uel (x —?t) I ” a a | - 
== 4 , dx 5 — x nm + >i x(x, t) (7.5b 
20 a 1 l — 
= m tSo = tS, ( (¢4.0C) 
2 t t 


Substituting So(1/t) and S\(1/t) from the appendix, we 
have 























(1 m)Cz,(t) = (o(t) + ¢t cos 06,(t) (7.6) 





= l l re 
Car, (t) es) = t(l — t) cos 8 ¢,(t) of t(o(t) (4.4) 
m 2 2 





where a = | corresponds to mid-chord. 


To obtain the response to an impulsive gust, we dif- 
ferentiate Eqs. (7.4c) and (7.5c) with respect to time, 








with the results 


O 7 10 C.,(t) = m[So(1/t) — (1/t)Ao(1/t)] (7.8) 














Cy,(t) = (m/2) faSo(1/t) — 4¢5,(1/t) + 
(2 — a) (1/t)Ao(1/t)] (7.9) 


Fic. 12 CLs Ut/c), the lift coefficient due to a unit impulsive 
gust, Vs. Tr. 
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Fic. 138. Cg! Ut/c), the pitching moment coefficient (about 
an axis two-thirds of the root chord from the apex) due to a unit 
impulsive gust, vs. 7 

(1/m)C,,(t) = cos 0(t) (7.10) 

i= l , g 

Cx,(t) pane ae — t}) cos 0 &,(t) + go(t) (7.11) 

\ O--2 4 6 8B 10 
The results [Eqs. (7.10) and (7.11)] are plotted in : , . P - 
Figs. 17 and 18. Fic 15. C( Ut/e), the rolling moment coefficient due to 

antisymmetrical gust of Fig. 14, vs. r 
(8) GENERAL DOWNWASH DISTRIBUTION; PITCHING 
AIRFOIL 

The general case of transient motion of an airfoil = 

may be specified by the downwash distribution 1O T 


D | 1(0)] se 
Wx, y,2) = — 2(x, vy, U)1(t) (S.la) 
° Dt P 
= —([z,(x, y, t) + a(x, y, OJ 1d — 
a(x, y, O)4(t)  (S.1b) 
where ¢ is the dimensionless time of Eq. (1.11). The 


response to the first two terms in Eq. (8.1b) can be 
obtained by Duhamel superposition of Eqs. (2.10) and 
(2.11), assuming 2(x, y, f) and 2,(x, y, t) to be continuous 
in both x and ¢, while the response to the last term can 
be obtained by calculating the response to —2(x, y, 0) 
l(t) with the aid of Eqs. (2.10) and (2.11) and then 
differentiating the result with respect to /, as in the 
case of the impulsive gust [see Section (5) ]. 

An example of some interest is the suddenly rotated 
airfoil for which 
(8.2) 


2(x, y, t) = ala — x) 
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Ci,( Ut/c), the rolling moment coefficient due to an 


antisymmetrical impulsive gust, vs. r 
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where a specifies the axis of rotation and a the angle. treated in Section (5). The response to the second term 
[he corresponding downwash distribution is given by may be obtained by substituting fp = 0 and w = (x — q) 
w(x, t) = all(t) + (x — a)d(t)] (8.3) in Eqs. (2.10) and (2.11), taking » = 0, and differen. 


The response to the first term in Eq. (8.3) has been _ tiating with respect to /—viz., 


fe) . x—&\. 
x5(x, t) = 2a tan pw tan 8 Pa Ay ; ) &(E— a) dé (8.4) 
0 


Carrying out the differentiation and then integrating by parts, we have 


2atanutand [{* x—é& : . P 
x3(x, t) = , : Ao ) 3& — 2x + avi + ax | dé (8.5) 


This must be added to Eq. (3.4) to obtain the overall response. While the result might be reasonably repre- 
sentative of a suddenly deflected control surface (of delta plan form), provided the time for deflection were small 
relative to (1 — sin #)~', we have not carried out any numerical calculations. However, the final results for the 
lift and moment coefficients can be expressed in terms of the S, integrals previously introduced. 


Appendix 
EVALUATION OF THE INTEGRALS S,(z) AND 7’,,(z) 


The required integrals are given by Eqs. (3.6) and (4.6) as 


5,(2) = / u"Ao(u) du, = 2S (Al) 
T(z) = / (1 — u)—"Ao(u) du, n = 2,3, 4 (A2) 
u2l(z—1) 

where Ao(u) is given by Eq. (2.6) and can be written 
Ao(u) = 1(2)Co(1 u) (A3) 
fo(1/e) = 1, 0<u<m, (mm = 1 — sin @) (Ada) 
= (1/m) cos~'! [(u — cos? 6)/u sin 8], UW <u < Ue (A4b) 
= @, u> us, (us = 1+ sin A) (Adc) 
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Fic. 17. Cr,(Ut/c), the two-dimensional lift coefficient due to Fic. 18. Cus( Ut/c), the two-dimensional (mid-chord) moment 


an impulsive gust, vs. 7. coefficient due to an impulsive gust, vs. r. 
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TRANSIENT LOADING OF WIDE DELTA AIRFOILS 


We consider first the S,,(s); substitution of Ao(w) in Eq. (A1) yields 
(n + 1)S,(z) = 3"7}, O<s< uy 
= [,(2), UW <3 < Ue 


= I ,,(u2), Z> Ue 


- ie 
T(z) = uy"t' + (n+ » f nro ( ) a 
J u 


Integrating /,(s) once by parts yields 


cos @ 


T,(2) = 3"**eo(1/z) + / u”"(— cos? @ + 2u — u*)” ? du 
ul 


T wo 


The last integral is elementary ;’ carrying out the evaluation and substituting in Eq. (A5), we obtain 
Sols) = 2f0(1/z) + cos 06,(1/2) 
2Si(z) = 37fo(1/2) + cos 06:(1/2z) — 2fe(1/2) 


3So(s) = 23fo(1/z) + [1 + (1/2) sin? 6] cos 6f,(1/2) — (1/2) (3 + 2)2f2(1/2) 


3 i. * 5 e 
£S3(3) = 2*fo ( ) + (1 + <sin? a) cos 06) ( ) — ( + ~ sin? @ + @ 2 + a) ooo ( ) 
\s 7 Z ) ) o Z 


oO 


(1(1/z) = 0, 0O< sz< mH 
= (1/7) cos! [(1 — 3z)/sin 6], Wy <3 < te 
= # Z> ue 
fo(1/z) = 0, GU <3 << ms 
cos @ ( “\(" )| n 
= — — | ‘ “<8 <8 
wT z Z 
= 0 ZD> Ue 


(n — 1)T,(z) = (1 — 3)°"" —1, O0<s<cuw 


QQ 


. l 
[,*(2) = (1 —m)""'-1+(am- vf (1 — un) "co ( Jaw, <3 
u1 u 


7 


, I 
-n-1 f (1 — uw) "vo ( ) du, ‘>i 
ar i u 


Integrating /,,*(z) once by parts yields 


_— I cos 6 
[,*(z) = (1 — 2) "fo —1(1 —s)+ 


ba us 


/ u—"(1 — u)'~"(— cos? 6 + 2u — u2)~ du 


0) 
>0) 


w 


u 
Introducing the partial fraction expansion 


n 1 
u—(1 —u)'""=u-'+ DS (1 —4)™ 
m=1 


(A5a) 
(A5b) 


(A5c) 


(A6) 


(A7) 


(AS) 
(AQ) 


(A10) 


(Al11) 


(Al2a) 
(A12b) 
(A12c) 


(Ala) 


(A13b) 


(Al3c) 


(Al4a) 
(Al4b) 


(Al4e) 


(Al5a) 


(A15b) 


(A16) 


(A17) 
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we have 


oe | cos 6 ” , : , 
I,*(s) = (1 — 3)" “bo * feo lil — 32s) + u-(— cos? 6 + 2u — u?) du + 


Tv l= sin @ 
n—1 "s 
cos 6 . . 
> / (1 — u)~"(—cos? 6 + 2u — u?) du (AIS 
T m=1 1+ sin @ 


where the upper and lower signs correspond to z less or greater than |, respectively. Carrying out the first inte 
gration’ and introducing the change of variable 


u= 1-—sin@cos¢ (419) 
in the remaining integrals, Eq. (AIS) reduces to 


" | cos 6 4 I “cos! ((1 3) sin @| : 
Tn*(2) = (U1 — a? = 1 }&o (:) + » (ese 6)” / (sec ¢)” dy (A20 


T m=1 1(s—1) 
Carrying out the last integrations’ and substituting the results in Eq. (A14), we obtain 


T2(z) = [2/(1 — 2)] [fo(1/2) —o&3(1/3)] (A21) 


) 
= ‘ me ii 5 me i zs i! 
273(2) = [( l al Z) 7 l l&o (‘) + | az cse¢* AC (‘) me | = o(-) 


Z 7 a 2 : I l 
= ; om + ese’ OF» = $s (A22 
7 l 1 2(3 — 23) Oa l | Z | 
37,(z) = [(1 — 2)~* — 1]fo ( -+- (1 2 esc? 0 ¢ ( — (1 + , ese’ 6 ca (A23) 
Zz 2 (1-3)? g 2 —s \s 
where 
¢3(1/s) = —[(1 — 3)/s] (cot 6/7) cosh™ (sin 6/1 — 2) (A24) 
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Comparison of Theoretical and Experimental 


Results for the Turbulent Boundary Laver 
in Supersonic Flow | 


* 


Along a Flat Plate 


CARLO FERRARI* 
Cornell Aeronautical Laboratory, Ine. 


SUMMARY 


On the basis of experimental results it is found necessary to 
restrict previously obtained theoretical data for the boundary 
layer arising from supersonic flow along a flat plate (placed 
parallel with the free-stream velocity vector) to an “‘inner’’ region, 
while a different type of solution is formulated for the ‘‘outer”’ 
region of the boundary layer, the characteristics of which depend 
ippreciably on the distance in the direction of the flow. The 
theory now conforms with the experimentally determined fact 
that the drag coefficient decreases with increasing Mach Number 


INTRODUCTION 


I A PREVIOUS REPORT,' which was the beginning of 
an investigation being directed toward a_ wide 
understanding of the properties of boundary layers in 
supersonic flow, the velocity distribution within a cross 
section of the boundary layer which arises from con- 
tact with a flat plate was determined for the case of 
turbulent flow; also, the resistance coefficient, C,, was 
obtained in terms of the Reynolds Number, R,. These 
deductions depend upon a parameter that had been 
considered as a universal constant and which was de- 
noted by 8. In the above-mentioned first report of 
this series, the presumed value of 6 was also given. 
This value was based on the assumption that for the 
same Reynolds Numbers, the friction coefficients, at 
least up to Mach Number = 2, would be equal or ap- 
proximately equal. Recent experiments by Young 
and Wilson have shown, on the contrary, that there 
is an appreciable decrease in C, with increasing 1/,. 
With use of these results as a starting point, the value of 
3 is now recalculated. Furthermore, from the same 
experimental evidence, it appears that the assumption 
upon which the preceding investigation was based 
Le., that the variables characterizing the physical and 
kinematical state of the fluid are independent of the 
\-coordinate along the plate—is far from being verified 
for the range of Reynolds Numbers within which the 
experiments were conducted. The investigation, there 
lore, is now extended by rejecting this restrictive as- 
sumption. To accomplish this purpose, a method is 
applied which is similar to that described in another 
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report’ of this series concerning the turbulent boundary 


layer. When this extension is carried out for the cases 
investigated here, the calculated results agree satis 
factorily with the experimental ones. 
SYMBOLS 
R = Reynolds Number of the undisturbed stream (R, = 
O.x/v 
/ = absolute temperature 
l = component of the velocity along the x-axis 
V = component of the velocity along the y-axis 
i = VU2+4 jy? 
p = density of the fluid 
M = viscosity coefficient 
v = kinematic viscosity 
p = pressure 
T = tangential friction stress 
rT, = tangential friction stress at the wall 
U* = friction velocity [U’* = (7,/p,) /?] 
, = limit velocity of the undisturbed stream when expanded 
intoa vacuum 
Ua = Ua/Vi 
f = U/U* 
f = U/l" 
0 = 0/0, =f/f 
a* = U*/V, 
= Mach Number 
y = T/T, 
6 = thickness of the boundary layet 
t = time 
( = coeflicient of local triction drag [c, = tp/(1/2) pala? | 
C, = coefficient of total friction drag for a plate of length x 
h = factor of proportionality relating friction stress to the 
product of the parameters p, 6, Oa, and 00 /dy 
Fs = x/L; n = y/L, x being taken in the direction of the 
wetted surface and where L = reference length 
Pp, = density in the external stream at a cross section of the 
layer having an abscissa equal to a 
= |see Eq. (9) of text | 
CU, = velocity of the external stream (its magnitude is taken 
to be equal to the component in the direction of the 
y-axis of the actual vector) 
Z = |see Eqs. (6’) and (6”) of text | 
y = ply = stream function (dimensions of a viscosity 
coefficient ) 
vi = ¥/paVol. 
6 = 6/L 
” o0\ y 
é* = 5 (: -—> )a = 
0 Pa Un 0 
g* = S*/L 
f 7 c C ) y 
6** = 6 = l1— = ])d 
0 P.U, U, 6 
a me Or 
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Bars and primes will be used to denote, respectively, the mean (1) While in the neighborhood of the wall the scale 
values of the above quantities and their fluctuations about their of the turbulence must essentially depend on the dis. 
mean values due to turbulence. The subscripts 0, p, and a cor- E : dts ie : 

: : : tance, y, from the wall on which the friction is oper- 
respond, respectively, to stagnation of the undisturbed stream, the : ee f 
wall, and the undisturbed stream. ative, in proximity to the external stream it seems 
likely that this scale is dependent on the thickness, 6, of 


the boundary layer. 


(A) DETERMINATION OF THE LAW OF VELOCITY 4 s ; ae 
: (2) For small values of y, the characteristic velocity 

DISTRIBUTION WITHIN THE BOUNDARY LAYER ; : : 
upon which the correlation among the turbulent fluc- 


tuations of the main velocity itself depends is U’,*, but 


for values of y close to 6, this characteristic velocity js 


Through reference to the first report! of this series on ct ley © ij 
§ ikely to be Uy. 


(A-1) Turbulent Friction Stress 


the turbulent boundary layer, it appears that the turbu- (3) In proximity to the external stream it is likely 
lent friction stress within the region of the boundary that r+ vanishes along with 0 /dy; on the other hand, 
layer close to the wetted wall may be written as for this same region, the variation of the several physi- 
ou U,* 7 cal quantities is certainly always related in some per- 

tT = pl a. es +8 — ppU,* U (1) ceptible way with the x-coordinate. At the wall, on 

4 : the contrary, at least starting from some suitable value 


of x, it seems likely that this dependence on the yx- 
coordinate is negligible in comparison to that related 
os ry ; to the y-coordinate. In this latter region close to the 

r = phiU, (OU /dy) (2) ee ; .- 
wall, therefore, the considerations still hold that lead 


while, for the region of the boundary close to the un- 
disturbed flow, it is true that 








one to assume that 8 is a constant for the “‘inside”’ of 
the boundary layer. However, within the ‘‘external’ 
region (which is external to the ‘‘inside’’ layer), this 


as was explained in the second report?’ of this investiga- 
tion of the turbulent boundary layer. 
The replacement of Eq. (2) by Eq. (1) may be justi- 


fied in the following way : assumption is no longer acceptable. 


A-2) Velocity in the ‘Inside’ of the Boundary Layer in Proximity to the Wall 
For sufficiently small values of y, 0p /Ox is equal to zero, and it is permissible to write 07/Oy = 0. Consequently, 
the distribution of velocity within this region of the boundary layer is still given by Eq. (22) of the first report 


of this series. The relationship is rewritten here. 


0 fi—Vyu*f V yurf t+ fe 
a EE log” : TR": log ust + fe 4 
c | (fi + fo) [Vv (u*/c) — fil Si (fi + fo) [V y(u*/o) + fe] Je 
l 
log (1 — ft = log y* + B (3) 


[V y(u*/c) + fe] [VW y(u*/c) — fil 


In the simpler case of a thermally insulated system, Eq. (3) may be reduced and becomes, upon reduction, 


OF 3 | mgt - f+ - : log (1 + fu*) + log (1 — cf), 
» oe og | _ Te og { ue | ee og | — Ci — 
cl. 2 [(u*/c)—-1] | ’ 2 [(u*/c) + 1] [(u*2/c2) — 1] wry 
log y¥* + B_ (3’ 
where c = Bu*. 





(A-3) Velocities Within the Boundary Layer in the Neighborhood of the External Stream 
The equation of motion for the boundary-layer flow in the neighborhood of the external stream is 
(0/dE) (aU2) + (0/0,) (BUT) = (0/0n)r 


where r is given by Eq. (2). Employing and y, as independent variables, one obtains 


aU (0/0t)U = (pU/p,U,) (Or O~1) (4 
On the other hand, it is true that 
O re) ee eS; fa) a ee EY a 
~ a ( a is — ) - ( andi ee “4 ) | 
, WwW On Uy oY Ba” UO Us 


Thus, it follows that 
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scale at 2 pUOo U 1 — u,7 0Z 5) 
i a -— = = (0 
de. Ba? UO Cg ta? OY 
)per- ; 
ems provided one sets therein ; 
6, of "U0 /0a =2 Uy U *o 
) ; : Dp , sat odo ; 
Z = 2Uu,’ | —— OC m= Bao tl — ay | tees (6) 
; 0 Oe Os o [1 — u,? yo? + (y — le)? 
city 
fluc- For a thermally insulated system (y = 1), one has simply 
but : 
Z = B/iy (6’ 
ty ie p/P ) 
For y ~ 1, one obtains 
ke . P = ” 
ety Z = (p/pa) + (¥ — 1)5(e) (6") 
and, 
1YSi- where 
per- 
. on - —] (y — 1)/2 — u,2 yo y—-1] Ua” ¥ 
alue ead’ RE 2 | 2.2 9 |+ 2 x 
ane Sia.* y+ Ky — 17 iy Ll — w,° yo" + Cy — De - Ua? ¥ + [(y — 1)?/4] 
© s » ° , 
I \ V ity y + [Cy — 1)2/4] — (vy — 1)/2 4+ 4.2 yo 
ated . log _ 
the 2V ua? ¥ + [(y — 1)2/4] | Vu y + [(y — 1)?/4] + (vy — 1)/2 — a? yo 
lead V ua? y¥ + [(y — 1)? 1] — (y —1)/2] —] y—1' (3 p 
” of log E = - = a - ) re (: — | - i. Fe ps + 
a V ua y + [ly — 12/4] + (vy — 1/25 2, ue? y + [(y—1)?/4]} ~ Pp Pp 
lal 
this te? 4 P V ua? y + [(y — 1)?/4] — (vy — 1)/2 + 4.’ v0 
~10¢g —_ 
945, 2 2/4]t? 5 : is . 
24 ta? y¥ + [(y—-1)7/4)5 7 | Va? y + [Cy — 1)2/4] + (vy — 1)/2 — 2? yo 
; V ua? y + [(y — 1)?/4] — (y — 19/2! 
og . 7 
V ua? y + [(y — 1)?/4] + (vy — 19/2) 
It results that 
tly, | ; 
ort! O7r/Oy, = pahd (Uq?/2u,g?) (1 — uy”) (07Z/Oy;7) 
On the other hand, it is true that 
U(0/dE)U = (U,2/2a”) (Ga?/p?) (1 — ta?) (0Z/0E) 
so that one has 
(9) B..pr0Z pO... dz 
02 = aS (8) 
Pa p o& Pa i OY; 
that is, 
OZ /de, = (p2/p,”) (U/Ug) (0°Z/0~*?) (8’) 
(3! ded 2 ae 
where y,* is set equal to ¥i/V hk and & = | 6dé. 
0 
Also, if one writes A-4) Solution of von Mises Eq. (9’ 
' 
ee Z (9) For the region of the boundary layer close to the ex 
(4 ternal stream, within which ‘“‘outer’’ region p > j, and 
where Zp = 1 + y(y — 1)8(1), then, U = U,, one has 
dz _ B? U O22 (9’) 0z,/0% = 072,/O~*? (9’’) 
oé i U,0 "e . ° ° — " 
™ m v1 together with the following boundary conditions: For 
which is formally analogous to Eq. (23) of the second y* = ©, 3, = 0; fory* = 0,2 = HT. Here, H de- 
report’ of this series. notes a function of & only, which will be specified 
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It will be presumed that it varies so slightly 
The solution 


shortly. 
that one may assume it to be constant. 
of Eq. (9”), satisfying the above conditions, is given by 


2 = Herfe (W,*/2V £,) (10) 


If it is now assumed that the “‘external’’ solution, Eq. 
(10), is a first approximation solution of Eq. (9’) itself, 
at least within a certain region, the exact solution of 
Eq. (9’) can be obtained by means of the method of 
successive approximations® in an analogous manner 
to the procedure suggested by von Karman’! in carrying 
out the study of the laminar boundary layer along a 
flat plate for the case of incompressible flow. 

Let the variable ¥,*/2V & be denoted now by means 
Also, note that Eq. (9’) can be 


written in the normal form, 


 d% ! | 1 Oz 
O€, oy,*? A(z) O£; 


where A(z) = (p?/p,”) (U/U,). 
In addition, let z, be the solution of Eq. (9’) which 


of the symbol &. 


corresponds to the (7 — 1)th approximation, so that 
z, is defined by the following equation: 


Oz, 0°, | | 1 O21 
O§; oy,*? A(Z,—1) O§, 





Remembering that it is true that A,z = 0 for ¢ = 0, one has C,’ = 0; 
¢ = o itisinferred that 
9 7 - 
\ Tw 0 wr 8 


Furthermore, one finds 


. 


; \ 
/ e ac f A, aGe' di, = = 
ee 0 Z 


and thus, 


Via , 
A,2 = erf ¢ A, 1Ge° 
2 0 
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J * 
0 
erfe ¢; dé, — / A, iGé ‘(ert ¢ — ert cas 11) 
0 
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If, in addition, one now sets A,s = 2, — 2,-1and A,_, = 


A(z,_1), then it follows that 


© Aes h (~ 1) Se 4 
of "(OW NA, OF 
( l 1) 9 
A,-2 Of, 
But, if z, is only a function of ¢, 2, = 2,(¢) (and be- 


cause it is true that 
02,/0& = —(¢/2&,) (dz,/de) 
02, /Oy~i*? = (1/4) (d2z,/d¢?) 
Cz. 1 Of; = —(¢ YE,) (dz, I dc) 
ete.), it may be inferred that 


d* d 
Anz + 2¢ 


a A,2 , A, iG = 
dc d¢ 


where 


/ A, sGe' dt, + C,' 
J 0 


> 


Ge’ dé) = / A, Ges eric dc, 
J 0 


Ge! (erf € — erf €))dé) 


This equation allows one to evaluate the nth approximation if the (7 — |)this known. 


The sequence of the functions 2; = 2, %, ete., cer- 
tainly converges,’ at least for values of ¢ confined within 
a suitable interval. Accordingly, the limit function of 
this sequence gives within this suitable interval the 
exact solution of Eq. (9’). The value of the function 
I], for each &, is to be determined by keeping in mind 
that the solution of Eq. (9’) which has just been ob- 
tained must produce velocity profiles for the series of 
cross sections along the boundary layer, which may be 


joined to the equivalent profiles afforded by Eq. (3) 
for the “‘inside’’ region of the boundary layer. One 
can easily see that a satisfactory “join” or connection 
and a satisfactory agreement with the experimental re- 
sults may be achieved by assuming for // values that, 
though changing from cross section to cross section 
along the boundary layer, vary so slowly as to justify 
the assumption made at the beginning of this Sec- 
tion. 


while from the equation A,z = 0 for 
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TURBULENT BOUNDARY L 
(B) LAW OF VARIATION OF THE THICKNESS OF THE 
BOUNDARY LAYER AND OF THE FRICTION COEFFICIENT 


B-1) Law of Variation of the Thickness of the Boundary 

Layer 

From the velocity profiles for the cross sections of 
the boundary layer, the thickness 6 of the layer can at 
once be inferred. It is worth pointing out that, even 
assuming the velocity distribution to be the one cor- 
responding to the external solution [Eq. (10)], one ob- 
tains, at least in the case of a thermally insulated sys- 
tem, a variation of 6 which is approximately linear with 
¢. Infact, itis true that 


BU, 
r= | 7 w= 2VEVEX 
J 8 


In addition, according to Eq. (6), one has 
B/B. = 1 —2 = [I — (C.2/V,))/(1 — (07/V,)) 


whence 


0/0, = V (ug? — 2)/[ue2(1 — 2)] 
ifu,” = U,?/ V7", and, consequently, 
(Ba/B) (Ue/U) = ta/V (1 — 8) (tg? — 2) 


Because of the assumption that // varies but slightly 
with £, one can consider the value ¢) of ¥i*/2V & for 
which z reaches a value corresponding, for instance, to 


(7 U', = 0.995 to be independent of & Consequently, 


: : dé 
6=2Véi, Vi u, = 
J0 V(1l — zs)(u,7 — 2) 


2V it, Vhu, BUI, u,) (13) 


wherein the following definition of B is employed: 


°$ dé 
BA, t,) = / (14) 
J0 V (1 — 2) (u," — 2 


Substituting in place of 


£, the integral form of this 
quantity, squaring Eq. (13), and differentiating (still 


assuming // = const.), one infers that 
db/dé = ys 

or (15) 
6 = 2u,"hB*t 


Because of the slowness of variation of B with &£, this 
equation demonstrates what has been stated above 
concerning the linearity of 6. 


(B-2) Friction Coefficient 


The value of the total friction coefficient is given by 
the relation 


AYES iN 
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TABLE | 
Velocity Distribution at Station 7 


U/U, y” y n 

0.6 73.40 0.1907 mm. 0. 0002384 
0.7 250.70 0.6483 mm 0. 0008106 
O.8 922.10 2.3850 mm 0. 0029810 
0.85 1,835.00 4.7460 mm. 0 .0059330 


B,U, 
ive 0 
| (1 “ ~)d “  6) 
E 0 Ua 2V& 


which, in the case of a thermally insulated system, be- 


comes 


. IVhvVé ti Ua” — 2 . 
Cc, = | — 2 dé (16) 
E 0 u,"{1 — 3) 


(C) NUMERICAL APPLICATION: COMPARISON 
BETWEEN THE THEORETICAL FORMULAS AND THE 
EXPERIMENTAL RESULTS BY YOUNG AND WILSON 


For the purpose of checking the formulas given in the 
preceding Sections, they have been applied to the de 
termination of the velocity distribution within the 
boundary layer along a flat plate under the same con- 
ditions as existed in the experiments performed by 
Young and Wilson.® The final results are presented 
below. The intermediate steps through which the 
value of the constants appearing in the formulas 
was obtained may be determined in a straightforward 
manner and are not reproduced here. 

At the 7th Station along the plate used in the ex- 


periments, where the abscissa x = 31.49 in. = 799.8 
mim., one has (assuming y = 1) J/, = 1.695; R, = 


17.7 K 10°: B = 0: c, = 0.0019. Then one finds 


u, = U,/V, = 0.6041 


U* 0, = 0.03868 
U* UP U,% 17.7 & 108 : 
es = 0.03868 ——_ (0.6351)! 
Vv, Ue te He 19.98 


| 
0.02586 K 107-7 em. 


hence, y = 0.02586 X 10~? y* em. 

The universal constant B is taken to be 2.2, since 
this is the value it assumes from the relation 
log (1 + by*) = log y* + B as y* becomes large and 
upon selection of the value 8.93 for 6 (see the first re 
port! of this series). From Eq. (3’) and with the aid 
of the above data, Table | may be constructed. 

The values in Table | are plotted on the graph of Fig. 
|; the experimental results by Young and Wilson are 
also given. It may be seen that, for C/U, < 0.82, the 
agreement between the theoretical and experimental 
results is fairly good, while the curve corresponding 


to Eq. (3 


, 


) departs more and more from the experi- 
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Fic. 2. Variation of the stream function upon traverse of the 


boundary layer, as related to a velocity parameter. 


mental one for l/l, exceeding the above-mentioned 
limit. In the following equations, 


y f? U ‘ f 1 — u,? U ; 
= an = 3 ee | 
* Fee U, i o 1 — u,2 (U2/U,) U4 


1 — nq? (02/02) 


let u,” = 0.3649. By use of the experimentally deter- 
mined values of l’/U,, Fig. 2 is plotted, giving 2 as a 
function of ¥,. Now because the experimental data 


for 6 equal to the value of y for which U/U, = 0.995 
yield the relation 10° 5 = 0.04 + 1.545 &, one obtains 


f* = 102 & = 0.04 & + 1.545 (#/2) (18) 


whence §* = 0.8123, for & = 1, and 2V é,* = 1.804. 
Consequently, it may be inferred that (upon selection 


of the value 2 for ¢ = 10y, Vb 2V ‘t)*, which implies 
the value 0.016 for y,) the following value for V h ap- 
plies: 


Vh = 0.16/1.804 (2) = 0.04435 


Then, assuming /7 = 0.24, one obtains Table 2 through 
use of Eq. (10). 

The values in Table 2 are also plotted in Fig. 2 as q 
graph of z versus ¥;. The agreement between this plot 
and the experimental curve appears to be satisfactory, 
It also joins up satisfactorily with the graph produced 
by the “‘internal’”’ solution, so that it seems unnecessary 
to push the approximation any further by following the 
method described in Section (A-4). 

By use of the experimental values of 2, as well as 
through aid of the calculated points on the graph of z 
versus ¥;—which together describe the solution for the 
“internal” and “‘external’’ equations of the boundary 
layer (the theoretical and experimental data curves are 
joined together as indicated in Fig. 2)—the friction 








coefficient has been evaluated from Eq. (16’). It is 
found that, for the experimental curve, C, = 0.0024 
and that, for the theoretical curve, C, = 0.00235. 

At Station 2, located on the plate so that the abscissa | 
x = 11.921n. = 302.8 mm., one has | 


é= 03785; MM, = 1.695: K, = 6.78 X 16 
(and still assuming y = 1) 
8=0; U,/V; = 0.6041 
If it is assumed that the local friction coefficient is | 
c, = 0.00211, it follows that 
U,*/U, = 0.04076; y = 0.02425 X 10-? y* cm. 
Through reference to Eq. (3’), Table 3 follows. 


TABLE 3 
Velocity Distribution at Station 2 





U/ ie y* y n 
0.6 52.9 0.128 mm, 0.0001605 
0.7 169.0 0.410 mm 0.000513 
0.8 581.5 1.411 mm, 0.001764 
0.85 t 117.0 2.711 mm. 0.003390 





The values in Table 3 of U/U, are plotted as a fune- 


tion of y on the graph of Fig. 3. The experimental re- 
sults are also given. Again it can be seen that the | 
agreement is good for U’/U, < 0.8 and that, for larger 


TABLE 2 
Data Required for Calculating Friction Coefficient at Station 7 


¢ 0 0.05 0.10 

v1 0 0.0004 0.0008 
S 0.2400 0.2265 0.2130 
c 0.35 0.40 0.50 

vy 0.0028 0.0032 0.0040 
z 0.1489 0.1372 0.1151 
4 1.00 1.20 1.40 

V1 0.0080 0.0096 0.0112 
z 0.03775 0.02152 0.01145 


0.15 0.20 0.25 0.30 
0.0012 0.0016 0.0020 0.0024 

0.1997 0.1865 0.1736 0.1611 

0.60 0.70 0.80 0.90 

0.0048 0.0056 0.0064 0.0072 

0.0951 0.0773 0.0619 0.0487 
1.60 1.80 2.00 { 
0.0128 0.0144 0.0160 

0.00568 0.00262 0.00112 
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TURBULENT BOUNDARY LAYER 


(’,, the departure of the theoretical from 


values of [ 
the experimental results progresses somewhat more 
rapidly than was true in the preceding case for Station 


(. 

Following the same procedure as used in treating the 
data for Station 7, the values of y and 2, based on the 
experimental data, have been calculated. Then, from 
Eq. (18), one has §* = 0.1258 and 2V ¢,* = 0.7094. 
Let £.9* be the value of £,* that corresponds to § = 1 
and set 

a [ek Jp 

W=h V &10*/é& (19) 
By use of this equation, into which the above-mentioned 
values are inserted, one obtains ¥; = 2.5421 y. This 
(z, ¥;)-curve is drawn in Fig. 4 to stand as a representa- 
tion of the experimental results. 

Now, taking /7 = 0.26 for the station under consid- 
eration, Eq. (10) yields Table 4. 

The graph of 2 versus ¥;, which depicts the values 
just calculated, is plotted on Fig. 4. Again, for this 
station, the agreement with experimental results, 
within the ‘‘external’’ region of the boundary layer, is 
satisfactory; the ‘‘join’’ with the previously evaluated 
“internal” solution is also good. 

The friction coefficients obtained from Eq. (16’) are 


now found to be 
(Cy)exp. = 0.0028; (C,)cate. = 0.0026 


A typical velocity profile has also been deter- 
mined by a reversal of these procedures and the use of 
presently established data. A hypothetical station 
along the plate has been selected for which J/, is 
still to be equal to 1.695 and for which R, = 10° X 
30; hence, € = 1.695. For this station, however, no 
experimental data are available. As a consequence 
of the investigation pursued in the first report! of this 
series, it was found that, for the asymptotic solution 
of the problem concerning the turbulent flow along 
a flat plate, c, and R, are related to each other by an 


equation of the type 
1/Vc, = A log (c,R,) + B 


If one accepts the validity of this relation and if use 
is made of the assumed values of c, for the two values of 
Reynolds Number (R, = 6.78 XK 108 and R, = 17.7 X 
10°) obtained previously, one is led to the conclusion 
in the present case that c, = 0.00176. Corresponding 
to this value of c, are the quantities 
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Table 5 results. 


TABLE 4 


Data Required for Computing Friction Coefficient at Station 2 


¢ 0 0.10 0.15 
V1 0.0004 0.0008 0.0012 
g 0.2453 0. 2307 0.2163 
¢ 0.40 0.50 0.60 
V1 0.0032 0. 0040 0.0048 
2 0.1486 0.1247 0.1030 
1.20 1.40 1.60 
y 0.0096 0.0112 0.0128 


0.02332 0.01241 0.00609 
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] 
0 
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TABLE 5 
Velocity Distribution at Station Where — = 1.695 
0/U. y* y n 

0.6 95.0 0.2554 mm 0.000319 
0.7 338.9 0.9104 mm, 0.001138 
O.8 1,311.0 3.5232 mm. 0.004405 
0.85 2,681.0 7.2019 mm. 0. 009004 
0.9 5,651.0 15.1830 mm. 0.018980 


The graph of l’/U, versus y is constructed on the 
basis of these calculations giving the curve depicted 


in Fig. 5. By a procedure analogous to that previously 

















used, the value of &* is determined to be &* = 2.287 
and 2V ¢,* = 3.024. It follows that 
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Fic. 6. Variation of the nondimensionalized and ‘“reduced”’ 
stream function as related to a velocity parameter. 
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Vi = hh V f10*/8:* = 1.804~1/3.024 = 0.5964y, 


From the above given values of ("/(", and y, the evaly. 
ation of Y, and z has been carried out, and the re a 
graph is plotted in Fig. 6. This curve of z vs. y, is 
drawn in asa solid line. 

The values for the solution of the 


have been obtained by taking 


“external” equa- 
tion, given in Table 
IT = 0.23. 

The graph that describes the relationship existing 
between these quantities also has been plotted on Fig, 
6 (represented by a long-dashed line); the two solu- 
“internal” and the ‘“‘external’’ ones 
The joining curve is drawn in as 


tions—1.e., the 
join together well. 
a dotted line. The value of the total friction coefticient 
was calculated and found to be C, = 0.00222. 

Finally, Fig. 7 presents the graph of C, 
It is derived from the 


versus R,, 
drawn on a logarithmic scale. 
three pairs of values of C, and R, which were given 
above. From this graph one infers, at least within the 
restrictions of the approximation being considered, 


that it may be stated that 


1/V C, = A®* log (R,C,) + Bo* (18’) 


and it is evident that the value A* = 1.697 may be 
deduced from those data. On the other hand, upon 
calculation of the local friction coefficient [by means of 
the equation, 1, Vc, A log (R,c,) Bo], one finds that 
the following approximate formulas hold: 
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Fic. 7. Plot of the total friction coefficient versus Reynolds 
Number for various Mach Numbers 


TABLE 6 
Data Required for Calculating Friction Coefficient at Station Where — = 1.695 
c 0 0.05 0.10 0.15 0.20 0.25 0.30 
V1 0 0.0004 0.0008 0.0012 0.0016 0.0020 0.0024 
Z 0.2300 ).2170 0.2041 0.1914 0.1788 0. 1664 0.1544 
c 0.35 0.40 0.50 0.60 0.70 0.80 
V1 0.0028 0.0032 0.0040 0.0048 0.0056 0. 0064 
z 0.1427 0.1315 0.1108 0.0911 0.0741 0.0593 
c 0.90 1.00 1.20 1.40 1.60 1.80 
v1 0.0072 0. 0080 0.0096 0.0112 0.0128 0.0144 
3 0.04671 0.03618 0.02063 0.01097 0. 00544 0.00251 
‘9 2.00 
0.0160 


v1 


0.00108 
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TURBULENT BOUNDARY LL 


"ise , 
= . (1 — 2AV C,) 
1+ 2AV¢, (19’) 


\2 
Il? 


A = A*V 14 2A*VC, = A*(1 + A*VC,) 


Consequently, from these relationships one obtains 
A = 1.823. 

These results permit one to check whether the values 
assumed for c, in the three cases considered above agree 
with the values of C, derived on the basis of the above 
assumptions. 

Introducing into Eq. (19’) the value of A just ob- 
tained, one finds: for C, = 0.0022, 


c, = 0.0022 (1 — 3.645 & 0.0472) 
0.00184 (instead of 0.00176) 


for C, = 0.0024, 


c, = 0.0024 (1 — 3.645 & 0.0490) 
0.00197 (instead of 0.0019) 


for C, = 0.00276, 


c, = 0.00276 (1 — 3.645 & 0.0525) = 
0.00223 (instead of 0.00211) 


The asymptotic equation, Eq. (55), deduced in the 
first report! of this series may now be tested. It ts 


repeated here as 
— Uy”) (1/F) 
where 


r a l i ) 4 
F = _ og (1 — Uy 
x1 — Bp) 


Ua 


B 
log (1 + u,) + log (1 — au.) 
2(1 + B) . l+ & 


It is interesting to note that, if 8 is assumed equal to 0 
and ay = 0.394, the asymptotic solution yields the 
value F = 0.5050 or, thus, A = 1.925. As far as the 
calculation of the thickness of the boundary layer is 
concerned, it may be easily determined through use of 
Eq. (15). Since the value of z which corresponds to 
C/U, = 0.995 is z = 0.00571, it may be seen from the 
tables that, in all of the three cases studied, this value 
implies that the magnitude of ¢» is close to (> = 0.16. 
The value of z which applies to Station 7 is first selected, 
and then making use of it in the evaluation of the co- 
efficient B of Eq. (15), one infers that u4,B = 1.979; 


hence, 
5 = 2 X 1.9792 & 0.001967E = 1.5404 &K 10~2(£) 


This proportionality factor is in good agreement with 


the previously given value. 


AYER 
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(D) APPLICATION OF THE METHOD JUST DESCRIBED 
TO THE PREDICTION ON A THEORETICAL BASIS OF THE 
VELOCITY PROFILE IN THE “‘INSIDE’’ OF THE BOUNDARY 
LAYER AND DETERMINATION OF THE RELATED DRAG 
COEFFICIENT 


For the practical application of the above formulas, 
not to the analysis of the experimental results but to 
the theoretical prediction of the velocity profiles within 
the boundary-layer cross sections, as well as to the pre- 
diction of the related friction drag coefficients, it is 
necessary to make use of the method of successive 
approximations. As first approximation values for c, 
one may take, for instance, those that are evolved, for 
given Mach Numbers and Reynolds Numbers, by use 
of the asymptotic solution, provided these results are 
diminished in the same ratio as in the cases considered 
in Section (C). 


Analogously, in order to predict the law of variation 
of the thickness of the boundary, one may take as the 
first approximation 6 = 0.332 Vc, &. 
formally analogous to that corresponding to the asymp- 
The coefficient 


This relation is 


totic solution, as well as to Eq. (15). 
of — has been deduced by taking into account the re- 
sults of the calculations performed in Section (C). It 
should be remarked that u,B does not depend appre- 
ciably upon u,. By following the method suggested 
and applied in Section (C), the assumptions just dis- 
cussed afford the means of computing the velocity pro- 
files for the various cross sections of the boundary 
layer. The construction of such profiles permits one, 
in turn, to determine the value of the coefficient C, 
and the thickness 6 of the boundary layer. It will then 
be possible to obtain a second approximation for c, and 
for 6; consequently, such an approximation will enable 
one to determine the corrections to be applied to the 
previously obtained results—that is, these corrections 
are to be taken as the differences between the assumed 


and the calculated coefficients. 


(E) GRAPH OF THE DRAG COEFFICIENTS DERIVABLE 
FROM THE ASYMPTOTIC SOLUTION 


The experiments by Young and Wilson clearly show 
a decrease in the drag coefficient with increasing Mach 
Number, and a comparable theoretical check on this 
coefficient may be obtained by use of the asymptotic 
equations of the flow. In addition, it appears from the 
preceding paragraphs that, in order to bring the theo- 
retical equations and the experimental results into 
satisfactory agreement (as far as the determination of 
the velocity in the region of the boundary layer close 
to the wall is concerned), it is necessary to set 8 = 0. 
With this provision in mind, a calculation has been 
carried through by means of the asymptotic equations 
[Eqs. (12) of the first report! of this series] in order 
to obtain the values of the coefficients A and By which 
define the asymptotic drag law, 

















564 JOURNAL OF THE AERONAUTICAL SCIENCES—AUGUST, 1951 
TABLE 7 7 
Calculations for Obtaining Drag Coefficients from the Asymptotic Solutions 
— F ~ — A 0.4342B ‘ 
M, B=0 B= —0.2 8B =0 B= —0.2 B=0 B= —0.2 : 
0 0.394 0.394 1.795 1.795 1.647 1.647 
0.5 0.400 0.392 1.809 1.849 1.638 1.657 B 
1.0 0.418 0.402 1.851 1.928 1.614 1.649 
1.5 0.445 0.420 1.915 2.027 1.580 1.628 pres 
2.0 0.476 0.444 1.995 2.138 1.540 1.599 
com{ 
The 
pin 
1/Vc, = A log (R,c,) + B — A log 2A 0.394 and 11.7, respectively. With the use of these 
values, the asymptotic formulas giving A and B (which, 
The calculation was carried out for several Mach for \/, = 0, become formally equal to the equations 
Numbers, as well as for two values of 8—i.e., for 8 = 0 suggested by — Karman as applicable in the case of 
an incompressible fluid) turn out to be, when \/, = 0, Tur 
and for 6 0.2. The values for the experimentally identically the same thing as the von Karman formula. T 
determined constants aj and } are taken to be equal to The results of these calculations are given in Table 7, Cc 
L 
Royé 
Marc 
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Turbulent Diffusion from a Source in 
Turbulent Shear Flow 


J. O. Hinze 
Royal Dutch/Shell Laboratory, Delft, Holland 


March 22, 1951 


I’ THIS JOURNAL, Dr. Corrsin'! explained the observed skewness 
in transverse temperature distributions close behind a fine, 











heated wire set perpendicular to both mean flow and direction 
of the greatest gradient by assuming that this transverse tem- 
perature distribution is simply equal to the probability density 
of the lateral velocity fluctuations. His hypothesis was that 
the skewness in this density is the result of the lateral gradient 
in lateral turbulence intensity. 

However, there is another aspect, not hitherto considered, in 
diffusion by continuous movements. It is not entirely correct 
to put the transverse distribution equal to the probability density 
of the lateral velocity fluctuations. It is better to assume this 
transverse distribution to be equal to the probability density 
\(y) of the lateral distances y covered by the fluid particles that 
have passed the source (provided the effect of molecular diffusion 
may be neglected ). 

Let be a transportable quantity; then, according to this as- 
sumption 
V(0) (1) 


I'(y) (oO) = V(y) 


According to the Lagrangian method for describing the path 
of a particle, the coordinates of the particles at time ¢ are 


y =u 
x=(U+ ud 
whence 
n = y/x =v/(U + u) (2) 


If we assume u/U < 1, Eq. (2) may be written 


v u u? ui 
l l U? U3 


If the average value of y for a large number of particles is 





taken, the result is 





If, in this expression the effect of the triple- and higher order 
correlations is negligibly small, the value of 7 must be propor- 
tional to the shearing stress in the flow at the point of the source, 


/ 
7 =~ —(uv/U?) = r/pU? (4) 
Since by definition 
i a n 
, Y(n) I'(y) 
n= ~- dn = ——~ » dn 
, S (0) ro) 


ae 


the average value of 7 is identical with the ordinate of the cen 
troid of the '-distribution surface. 

The skewed [-distribution is determined by the probability 
density Y(y). Let F(u, v) be the probability density of the 
joint distribution u, v. According to the relation [Eq. (2)] be- 
tween y and the turbulence velocities u and v, the probability 
density Y(y) apparently is obtained by taking from F(x, v) all 
those values for which the relation [Eq. (2)] holds, or we must 


integrate F(u,v) along a pathv/(U + u) =» =constant. Thus, 
le 
YV(y) = F(u, v) ds 
along the pathv/(U + u) =n 
Since 
ds = V1 + (dv/du)? du = V1 + n? du 

the expression for Y(n) becomes 

Yin) = © hohe Flu, n( U + u)] du (5) 


The function F(u, v) has to satisfy the conditions 


1 Joe 
i= 7 F(u,v) udu dv = 0 


; bias 
= F(u,v)vdudv = 0 


¢ 
~. 


l . 
i. 2s 
> / F(u, v) u? du dv = u? 
U2 . 
I a 
> / F(u, v) v? du dv = v? 
U2 
i. ee 

F(u, v) du dv = 1 

U? = 


At the present juncture, nothing is known about the function 
But a function that satisfies all the above conditions 


F(u, v). 
is the normal correlation function, 
u ' 
F(u,v) = x 
2ru'v’ V1 — R?* 
—1 nu? 2R uv v2 
exp : _ “Erie. oer (6) 
21 — R*)\u" ut vy’? 
where 


R = i0/u'v’, uv =~ Vu, v = Vv’? 


If we assume the turbulence velocities to be normally corre- 
lated, though this implies that the marginal distributions or the 
probability densities of u and v separately are Gaussian, the ex- 
pression for (7) becomes 

’ 1 + 7? 
V(m) = [ om 
2r(v"? 


— 2Rn u'v’ + n?u"?) 


exp 


' | 1 + 9" 
I'(n) | rs a 
ro) 


if at the same time 


; cor 
obo 
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re . can 
T= —puv = —pRu'v the expectations that prompted these supplementary te: ; 
ate 3 t the : 
is introduced in the expression for I’. the new drag coefficients were smaller than those presented j: it th 
It may be interesting to see to what expression (7) reduces a L. Ms might pu Nae that the interference initiate 
when the turbulence is isotropic. In this case, r = 0 and in the flow a en the sphere by a rear sting support is at leas; 
Age a equally as serious as that caused by the cross-stream rod Tr 
date no optimum support technique has been developed here Tt 
i ' = exn | — [/2n2/ [2y’2 d 2)]t ° 1e 
P(9)/TP(O) = exp | — U%n?/[2v°(1 + 9?) J} (8) Also apparent in the more recent drag results are certain djs 
According to this result, the spread for great values of 7 is ¢repancies between drag-coefficient values obtained at a give 
greater than according to the Gaussian error law. This differ- Reynolds Number, but with different size spheres. These djs whet 
ence is caused by the fact that the effect of the turbulence ve crepancies were identifiable because of an improvement in pres as 
locity u on the spread has not been neglected as it had been SUre-measuring instrumentation which reduced experiment 
hitherto. scatter considerably. They have not been explained satisfa 
It is intended to present the results of applying the above con torily in terms of the small variations of Mach Number involve there 
siderations to measurements of the local diffusion of heat from a but may be due to these, in combination with such factors 
heated wire located in the turbulent mixing zone of a plane free — the nonuniformity and limited extent of the free jet, nonsimilarity 
jet to the general discussion on heat transfer and related fluid of support geometries, and variations of sphere surface conditioy Putt 
mechanics at the meeting of the Institution of Mechanical Engi- 
neers and the American Society of Mechanical Engineers to be REFERENCE 
Idi septe . 95 R Subs 
held in London, Septembe r, 1951. 1 Kane, E. D., Sphere Drag Data at Supersonic Speeds and Low Reynolj 
Vumbers, Journal of the Aeronautical Sciences, Vol. 18, No. 4, pp. 25 
REFERENCE April, 1951 
TI 
1 Corrsin, S., Hypothesis for Skewness of the Probability Density of the 7 
Lateral Velocity Fluctuations in Turbulent Shear Flow, Journal of the Aero = 
nautical Sciences, Vol. 17, No. 7, p. 396, July, 1950 , 
. . wher 
Effective End Restraint of Columns by 
- Frequency Measurements 
Putt 
Harold Lurie 
Note on Sphere Drag Data The RAND Corporation 
April 24, 1951 As tl 
Frederick S. Sherman Eqs 
University of California at Berkeley ¢ A RECENT INVESTIGATION, I came across a paper by Stephens 
May 4, 1951 in this Journal! which is based on a fallacious argument, anc 
apparently his erroneous analysis has not been pointed out sub 
N HIS RECENT REPORT! on sphere drag data at supersonic S¢quently. Stephens proposed that the end restraint of a bean B 
I speeds and low Reynolds Numbers, Kane remarks that the could be found by measuring the natural vibration frequency of ; 
experimental drag coeflicients obtained may have been fictitiously the beam, where the end restraint is defined as the coefficient 
low because of interference between the cross-stream support rod = 1! the basic Euler buckling equation Th 
_— ° . ° ° er 
and the sphere boundary layer. To investigate this matter, his Per = cr®EI/l? 
tests were repeated (again under the sponsorship of O.N.R. ae on I 
, 2 aoa lo correlate the degree of end fixity with the natural frequency 
O.A.R.) last summer, using the same wind-tunnel facilities but ‘ : : ae 
: : Stephens used the value of c for the following four different types 
employing a rear sting support arrangement for the sphere 2 . ' 
ad Pecks : “a of end restraint: (a) one end fixed, one end free (¢c = '/4); (b 
models. The principal results of these tests are shown in Fig. 1, —_ wor 1 va 4 | fixed 1 pi ' Hen 
é : ‘ 4 oe . 0th ends pinnec = 1); (c) one end fixed, one end pinned (c = | 
in which Kane’s results are indicated by the curve that best fits a she : tee ‘ site aie end f 
: 2.047); and (d) both ends fixed (¢ = 4). 
his data. The f as j , f i . of 14 
_ — P . ‘ le frequency equations for these four type restraint wer = 
The drag coefficients obtained with the rear sting support | ‘dit. = scab “dsl om : : fixity 
‘ . a then calculated, and the first eigenvalue determined for eae! . 
method differ strongly from those obtained with,a cross-stream Tt et : ' Hed tt ‘ concl 
° : “ase 1 square of the eigenvalue wi “alle le ‘‘Trequency > 
support only for Reynolds Number <100. Below this value “** “i — ta + Se oe Rs _ fiable 
: : “ constant,’’ denoted by A. The values for A are 3.516, 9.87 
the discrepancy becomes large, reaching about 20 per cent 15.421 1 29.37 hale i ; ' f B In 
‘ 5.421, and 22.373, respectively, for the above ur cases. By 
at the lowest Reynolds Numbers encountered Contrary to ee ere sis - : y ; sacar cl iat a $3 thrus 
plotting A against c, Stephens obtained a single curve. Thus, t 
determine the degree of end restraint, he stated that it was onl 
necessary to determine the natural vibration frequency by test wher 
after which A could be determined from the equation m fot 
SPHERES - TAILSTING MOUNTED f = (K/2m i?) WEI/p actu 
22<" < 2.6 sa ti stant 
! < Re < 1100 The corresponding value of c, which is the measure of the e! depe 
LEGEND: fixity, could then be found from the curve of A against c. abov 
0.10 IN. DIA. SPHERE However, his proposed relationship between A and c¢ is 1 ods, 
sy 4 . ‘i unique, since there may be more than one value of A correspond finall 
1.00 « ing to each value of c. This might be suspected from the physi 
1.25 TA of the problem. It is feasible for two identical bars with differ 
ently restrained ends to have the same resultant fixity coefficiet a 
c. However, different modes of vibration and, hence, differen As pe 
frequencies would be anticipated nautic 
A simple example will illustrate this. Consider a strut elas * ha 
. ‘ . Ae 
tically restrained at each end, and suppose the restraint is suc pes 
: i 56 -_s ake 
10 100 1000 as to give an end fixity coefficient equal to 2.047 (the same asf0) Comp 
REYNOLDS NUMBER FIG. | one end fixed, one end pinned). The value of K for such a bea 
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Suppose the spring constant at each end is 
Taking the origin of coordinates 


can be calculated. 
the same, and denote it by a. 
it the center of a beam of length 2/, the boundary conditions are 


v= 0 atx = tL | a 
) 
Elv” = av’ atz = tL } 
The general solution of the buckling problem is? 


= A cos Bx + B sin Bx + Cx + D 


where 6B? = P EI. 


As 
P = crEl/AL? 


therefore, 


bo 


B2w-arVy c/2L & 
Putting in the boundary conditions |Eqs. (1) ], 
(EI/a)B = tan BL 
Substituting the value for 6 from Eq. (2), 
a = (wEIV c/2L) cot (r¥ c/2) (3) 
The general solution for free vibrations is* 
= A cos \¥x + B sin Ax + C cosh Ax + D sinh Ax 


where 
? = wy p El 

Putting in the boundary conditions |{Eqs. (1) |], 

a = 2 El)d/(tan AL + tanh AL) (4) 
As the value of a is to be the same for buckling and vibrations, 
Eqs. (3) and (4) may be equated Therefore, 

2rAL 

tan AL + tanh AL 


cot 


wrV ¢ wrV « 
5 3 


By taking c = 2.047 as for the fixed-pinned strut, Eq. (5) becomes 


2,\L/(tan XL + tanh AL) = —1.8052 


Therefore, 
2\L = 3.776 
K = (2AL)? = 14.258 
Hence, for a symmetric strut with elastically restrained ends, an 
end fixity coefficient of 2.047 corresponds to a frequency constant 
of 14.258; whereas for a fixed-pinned strut with the same end 
fixity coefficient, the frequency constant is 15.421. It may be 
concluded that Stephen’s single curve relating ¢ to A is not justi 
fiable. 
In the same paper, Stephens derives the relation between end 


thrust and vibration frequency as 
P =m Per {1 — (f/fo)*] 


where m is a constant depending on the end restraint. He lists 


m for various end conditions. This again is incorrect, since m is 
actually a function of P if the ends are not pinned 
stant m was determined by the energy method, which makes it 


According to the 


The con 
dependent upon the approximation used. 
above formula, P = Per when f = 0. 
ods, it can be shown that as P is increased, f decreases until it 


However, using exact meth 


finally vanishes at P = Per. 
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Comment on the Calculation of Deflections of 
Indeterminate Structures 


Carlos Thalau 
Instituto Aerotecnico, Cordoba, Argentina 
May 8, 1951 


wr REFERENCE TO Note on the Calculation of Deflections of 
Indeterminate Structures by Bruno A. Boley and Roland 
H. Moore, published in the JoURNAL OF AERONAUTICAL 
ScrencEs, August, 1950, and the considerations of H. F. Michiel 
sen in the same publication, April, 1951, concerning the same 


THE 


item, the following remarks may be of interest 


The fact that in calculating a deflection 6;%’ of an indeter 


minate structure it is possible to write 


i.e., to substitute for the indeterminate structure a conveniently 
chosen determinate structure in one of the two loading systems by 
applying the principle of virtual displacements (prime for the 
indeterminate structure and sub zero for the determinate struc 
ture) has been known as the “law of reduction” (Reduktions 
satz) fora number of years 

As I remember, Mueller-Breslau, in his Graphische Statik der 
Baukonstruktionen (published by Kroener, Leipzig), and Kammer, 
in his Statisch unbestimmte Hauptsysteme (published by Springer, 
Berlin), have already treated this fact in detail. Still other 
sources probably exist 
Flugzeugstatik [Thalau-Teich 
mann (published by Springer, Berlin, 1933) |, which I had the 


In our book Aufgaben aus der 


honor to present to your Institute at that time, we treated the 
law of reduction with a few examples on pages 202 and 285, since 
it belongs among the tools of every engineer concerned with 
statics of indeterminate structures 

In agreement with the wishes of the authors of the note, the 
“law of reduction” already has found its way into a few books 


Note on ‘‘Heat Transfer by Laminar Flow 
from a Rotating Plate’’* 


Arthur N. Tifford 
The Ohio State University, Columbus 
May 8, 1951 


— RECENT ABSTRACT! of an exact solution of the energy 
equation in,the case of laminar flow along an isothermal 


rotating plate is interesting. It is related to an exact solution 


of the energy equation in the case of flow near a forward stagna- 
tion point obtained here last summer in work supported by the 
Office of Air Research, U.S.A.F. The latter analysis resulted 
from a desire to investigate the influence upon boundary-layer 


solutions of extremely large surface temperature variations 


under such conditions, Prandtl’s order-of-magnitude 


arguments do not significantly simplify the equations, the exact 


Since, 
equations were analyzed. The detailed results will be reported 
soon. 

Our earlier exact solution suggests that the rotating plate being 
considered need not necessarily be isothermal in order for the 
complete energy equation to be tractable. A little thought on the 
matter reveals that the form of solution given in reference | may, 
as to be applicable 
his is 


in fact, be interpreted in such a manner 
when the surface temperature distribution is parabolic 
only a special case, however, of the more general surface tempera 
ture distribution 


* Work supported by the Office of Air Research, U.S.A.F 
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. don?” 


Tw = Tret. + do + aor? + art +.. 


that may be subjected to an exact analysis. The associated 


general solution may be written as, 
y = T ambient = G(s) + Aogo( z) + d2go() + asqa(z) + 
. Aongon(Z) + 1? [5(s) + a2s(s) + aasi(z) + 
‘ AonSon(2) | + r*[asts(z) + coe Aonton (z)] + 
ooo T 2" aon yon(2) 


where all functions are simply related to universal functions (for 


CAL SCIENCES—AUCUST, 


a given Prandtl Number) and where the barred functions com. 
pletely account for the contribution of frictional heating, ]f 
would be of additional interest for Drs. Millsaps and Pohlhausen 
to consider some of the facetst of the more general solution ig 


their forthcoming paper. 

T For example, as in the case of an isothermal spinning surface,’ there js 
an analogy, when the Prandtl Number is unity, between the rate of heat 
transfer and the power expenditure required to overcome frictional resistance 
to rotation In the case of the rotating plate, however, the analogy applies 
even when the surface temperature distribution is parabolic. 
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Similarity Laws for Slender Bodies of Revolution 
in Hypersonic Flows 


(Continued from page 330) 


similarity law for three-dimensional flow has been 
presented, including its application to finite span wings 
and wing-body combinations. However, the accuracy 
and range of applicability has not as yet been es- 
tablished. It seems clear from these results that future 
effort should be directed toward further verification of 
the law in its various forms and toward determination 


of its range of applicability. 
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